AD-A104  933 

UNCLASSIFIED 


SYSTEMS  TECHNOLOGY  INC  HAWTHORNE  CA  F/G  _/0  . 

USER'S  OUIDE  TO  BIOOYN-BOl  AN  INTERACTIVE  SOFTWARE  PACKAGE  FOR  — ETC  KJ 
JUL  81  S  A  RIEDEL*  R  E  MAGDALENO*  H  R  JFX  F3361S-79-C-051Q 

STI-TR-1146-1  AFAMRL-TR-01-59  NL 


AFAMRL-TR-81-59 


i  AD  A1  0493^ 

V 

l  '* 


USER’S  GUIDE  TO  BIODYN-80; 

AN  INTERACTIVE  SOFTWARE  PACKAGE  FOR 
MODELING  BIODYNAMIC  FEEDTHROUGH  TO 
A  PILOTS  HANDS,  HEAD,  AND  EYES 


NOTICES 


/ 


When  US  Government  drawings,  specifications,  or  other  data  are  used  tor  any  purpose  other  than  a  definitely  related 
Government  procurement  operation,  the  Government  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever, 
and  toe  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way  supplied  the  said  drawings,  specifi¬ 
cations,  or  other  data,  is  not  to  be  regarded  by  implication  or  otherwise,  as  in  any  manner  licensing  the  holder  or 
any  other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented 
invention  that  may  in  any  way  be  related  thereto. 

Please  do  not  request  copies  of  this  report  from  Air  Force  Aerospace  Medical  Research  Laboratory.  Additional 
copies  may  be  purchased  from: 


National  Technical  Information  Service 
S285  Port  Royal  Road 
Springfield,  Virginia  22161 

Federal  Government  agencies  and  their  contractors  registered  with  Defense  Documentation  Center  should  direct 
requests  for  copies  of  this  report  to: 


Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 


TECHNICAL  REVIEW  AND  APPROVAL 

AFAMRL-TR- 81-59 


this  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (PA)  and  is  releasable  to  the  National  Technical 
Information  Service  (NTTS).  At  NTIS,  it  will  be  available  to  toe  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


FOR  THE  COMMANDER 


HENNING  E.  VON  GIERKE 
Director 

Biodynamics  and  Bioengineering  Division 

Air  Force  Aerospace  Medical  Research  Laboratory 


'.jir  - 


SECURITY  CLASSIFICATION  Of  THIS  PACE  (Whtn  Dm(m  Bnffd) 


j  |  /  REPORT  DOCUMENTATION 

PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

T 

RE^RTVjiUMBER 

AFAMRl|{TR-81-59  (  . 

2  GOVT  ACCESSION  NO. 

hb-/)J.cv(/ 

S  RECIPIENT'S  CATALOG  NUMBER 

S3 

4 

T*T fend  SUBffPrtJ  — -  ~~ . - 

USER’S  £UIDE  TO  ^0DYN-80;  INTERACTIVE  \  n 

Software  package  per^DELING  ^iodynamic  \  v  7 

FEEDTHROUGH  TO  A  PILOT'S  HANDS,  HEAD,  AND  ) 

1YES*  '  "  sZZI<t 

5  TYPE^B  -RBReRT  »  PeRtOO  COVERED 

Final  f(ep6rt*  J  u 

177 9  io  11/80 

J.  PERFORMING  ORG.  REPORT  NUMBER 

r-TR-1 146-1  / 

nuTHORr*)  ,  ^  /  /y 1 

: Susan  A.  Riedel 

Raymond  E.  Magdaleno  /  ■ 

Henrv  R.  Jex  j 

-p-  BORTRACT  «WGRANT  NUMBS Rf«) 

F33615-79-C-0519 

9.  PERFORMING  ORGANIZATION  NAME  ANO  ADORESS 

Systems  Technology,  Inc. 

13766  S.  Hawthorne  Blvd.  / 

Hawthorne  CA  90250  .  .  . 

10.  PROGRAM  el ewtnt.  project.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

61102F  _  .  j  - 

T2712I-V3-20 

It  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Air  Force  Aerospace  Medical  Research  Laboratory,  t 
Aerospace  Medical  Division,  Air  Force  Systems 
Command,  Wright-Patterson  AFB  OH  45433 

DATE 

'/)  July  1981 

»S.  NUMBER  or  PAGES  ✓  J 

73  y  ,  /  - 

14  MONITORING  AGENCY  NAME  ft  ADDRESS  (it  different  from  Controlling  Office) 

IS.  SECURITY  CLASS.  (oTlKTt  rtport) 

Unclassified 

15*.  DECLASSIFICATION/ DOWNGRADING 

schedule 

is  Distribution  statemen t  (ot  thlt  Rtpon) 


Approved  for  public  release;  distribution  unlimited. 


J  ’7.  DISTRIBUTION  STATEMENT  (of  the  mbatract  entered  in  Block  20,  It  different  from  Report) 


18.  supplementary  notes 


19  KEY  WOPOS  ( Continue  on  referee  aide  It  neceaamry  and  Identify  by  block  number) 

Biodynamics  Control  Performance 

Whole-body  Vibration  Man/Machine  Interface 

Visual  Performance  Computer  Program:  BI0DYN-80 

Human  Transmissabil ity 


JO  ABSTRACT  f Continue  on  reverae  aide  If  neceaaary  and  Identify  by  block  number ) 

This  is  the  User's  Manual  for* -BI0DYN-80^  an  interactive  computer  program 
for  solving  the  dynamic  motion  response  of  seated  aircrew  (pilots  or  equipment 
operators)  working  in  a  biodynamic  environment  (such  as  vibration  or  changing 
levels  of  acceleration  in  aircraft,  surface  vehicles,  or  motion  simulators). 
The  multi-degree-of-freedom,  lumped-parameter  model  includes  elements  for: 
pelvis,  torso,  neck,  head,  eyes,  upper  and  lower  arms,  hand-grip  and  control 
stick  (center-stick  with  arm  rest).  An  active  21-parameter  neuromuscular 
system  is  operable  with  either  the  limb  or  head  control.  The  eve's  image 

oo ,  1473.  EDITION  OF  I  NOV  «S  1$  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whtn  Dtlt  BnIPrtd) 


SECURITY  CLASSIFICATION  OF  THIS  PAOE(TTh«n  Palm  B nlmnd) 


20.  ABSTRACT  (continued) 

motion  dynamics  include  a  13-parameter  model  for  the  fixation  and  vestibulo- 
ocular  reflexes.  The  nonlinear  equations  of  motion  of  the  body's  main  elements 
are  linearized  about  a  general  set  of  body-limb-head  postures,  ranging  from 
prone,  to  erect,  to  supine.  Inputs  are  vertical  or  fore-aft  accelerations  of 
the  seat,  while  a  wide  variety  of  outputs  are  available,  including  motions  of 
the  shoulder,  head,  eyes,  arm,  hand,  or  control  stick.  Typical  parameter  ranges 
and  sources  are  given,  along  with  two  '‘typical*  sets:  a  seated-pilot  with 
center-stick,  and  a  seated  crewman  with  hands  in  lap. 

The  program  can  be  run  in  a  mode  interactive  with  the  user,  or  in  a  batch 
mode.  The  program  helps  the  user  load  the  necessary  parameters,  and  offers  a 
'•quick-look*  printer  plot  of  the  resulting  frequency  responses,  in  standard 
Bode-plot  forms.  The  program  also  can  produce  biomechanical  transmissability 
files  needed  for  use  in  the  Air  Force's  PIVIB  program  for  tracking  performance 
estimation. 

Procedures  and  examples  are  given  for  both  a  Cyber  175  version  available 
on  "INTERCOM"  for  Wright-Patterson  Air  Force  Base  users, -and  a  PDP-10  version, 
available  on  the  Tymshare  national  computer  network  for  other  users.  < 

y"  . 

I 


SECURITY  CLASSIFICATION  OF  THIS  PAGEftniwt  D*f  Entfd) 


PREFACE 


This  report  was  prepared  by  Systems  Technology,  Inc.,  Hawthorne, 
California,  under  Contract  F33615-79-C-0519,  for  the  Air  Force  Aerospace 
Medical  Research  Laboratory  (AFAMRL).  Mr.  Charles  Harmon,  Biodynamic  Effects 
Branch  of  the  Biodynamics  and  Bioengineering  Division,  served  as  the  technical 
monitor  for  AFAMRL  in  support  of  Project  2312-V3-20,  "Man-Machine  Interface 
Model." 


Accesric n  For 

"fits  gfm-t 


i  . -- 


.  A  Vf> 1  l 


i  v  f!o- 
;  .  :/cr 


TABLE  OP  CONTENTS 


INTRODUCTION .  4 

A  Scope . . . . ...  4 

B  Background .  4 

C  Applications . 5 

SOFTWARE  OVERVIEW .  6 

MODEL  DESCRIPTION .  10 

A  Biomechanical  Model . 10 

B  Limb  Neuromuscular  Model . 12 

C  Neck/Head/Eye  Models . 12 

USING  BI0DYN-80 .  17 

A  CREATE .  17 

B  BIODYN .  27 

C  Plot .  28 

D  Interface  with  PIVIB .  28 

E  For  Assistance... .  35 

APPENDIX  A  DEFINITION  OF  PARAMETERS,  SYMBOLS,  RANGES 

AND  TYPICAL  VALUES .  36 

APPENDIX  B  DEFINITION  OF  FORCING  FUNCTION  AND 

RESPONSE  VARIABLES .  45 

APPENDIX  C  EQUATIONS  OF  MOTION .  49 

APPENDIX  D  TYMSHARE  EXAMPLE  PROBLEM .  57 

APPENDIX  E  INTERCOM  EXAMPLE  PROBLEM .  66 

REFERENCES .  69 


2 


LIST  OF  ILLUSTRATIONS 


Figure  Pfl8e 

1  Overview  of  BI0DYN-80 .  7 

2  BI0DYN-80  Flow  Chart .  8 

3  Main  Biomechanical  Elements . .  11 

4  Linearized  Limb  Neuromuscular  System .  13 

5  Simplified  Head  Pointing  Neuromuscular  System .  14 

6  Block  Diagram  for  Head,  Eye,  Vestibular  Models 

and  Display  Kinematics . 15 

7  Example  Long  Format  File  Listing . . .  20 

8  Example  Short  Format  PARAMETER  File  Listing .  22 

9  Format  for  BI0DYN-80  Transfer  Function  Printouts .  29 

10  Format  Used  for  BI0DYN-80  TAPE 19  File  of 

Transfer  Functions .  30 

11  Example  Output  from  Plot . 31 

12  BIODYN-80  Transfer  Functions  Which  May  Be  Used 

as  PIVIB  Input*.. . 34 

LIST  OF  TABUS 

Table 

1  Parameter  and  Symbol  Definitions,  Ranges,  and 

Typical  Values . 36 

2  Parameter  File  Structure  (Short  Format) .  43 

3  Input  (Forcing  Function)  Variables .  45 

4  Output  (Response)  Variables .  46 


3 


ihtbodoctioh 


A.  soon 

BIODYN-80  is  a  versatile  computational  tool  used  to  determine  trans- 
missibilities  (transfer  functions)  between  vertical  and/or  fore-aft 
vibrational  inputs  and  important  biodynamic  outputs,  such  as  motions  of 
the  torso,  head,  eyes,  arms  or  hands.  The  program  scenario  assumes  a 
seated  pilot,  gripping  an  arbitrary- angle  stick  and  viewing  a  display, 
possibly  engaged  in  a  tracking  task.  The  physical  model  uses  an  "iso¬ 
morphic,"  (lumped  parameter,)  approach  to  represent  the  relevant  por¬ 
tions  of  the  whole-body  torso,  limbs  and  head,  as  well  as  postural 
compliances  among  the  joints.  The  implementation  of  this  model  includes 
a  chain  of  Interacting  parallel  and  serial  second-order  elements,  with 
neuromuscular  and  other  force  feedbacks  at  the  arm  or  head.  The  result¬ 
ing  equations  are  in  "second-order  element"  matrix  form  and  apply  to  a 
wide  range  of  seated  postures.  A  separate  input  file,  which  describes 
the  particular  set  of  parameters  to  be  used,  is  created  by  the  user 
(usually  by  modifying  one  of  a  cataloged  set).  This  file  is  incorpor¬ 
ated  in  the  matrix  to  produce  linearized  coefficients  for  perturbations 
about  the  selected  equilibrium  posture.  A  variety  of  outputs  and  inputs 
can  be  specified  to  evaluate  the  desired  transmissibility  transfer  func¬ 
tions,  and  these  are  written  to  a  file  in  formatted  form  for  plotting  or 
use  in  other  programs.  On-line,  printer-drawn  frequency  response  plots 
(Bode  format)  are  available  to  aid  in  interactive  user-computer  opera¬ 
tions,  or  to  screen  significant  results  from  batch  runs. 

■  .  BAGKOOUaD 

BI0DYN-80  is  one  result  of  a  several  year  small-scale  development 
effort,  reported  in  detail  in  References  1-6.  It  is  based  on  vibration 


Measurements  made  at  the  Aeroapace  Medical  Research  Lab/Blodynamics 
Division  (AMRL/BB)  and  elsewhere.  Most  of  the  torso,  limb  and  stick 
model  elements  are  based  on  independent  vibration  measurements  (e.g., 
Reference  5),  and  the  neck,  head  and  eye  effects  show  promising  corre¬ 
lations  with  the  few  available  measurements  on  image  motion  effects 
(Reference  6).  However,  many  aspects  remain  to  be  explored,  validated 
or  upgraded  as  more  experiments  are  run  and  interpreted  via  BI0DYN-80. 

C.  APPLICATIONS 

The  possible  applications  of  this  program  are  many.  It  should  be 
used  in  the  early  stages  of  experimental  design  for  determining  the 
optimal  locations  for  vibration  measurements  and/or  selection  of  fre¬ 
quencies.  It  can  also  be  used  by  development  engineers  for  solving 
practical  pilot-vehicle  interface  design  problems  such  as  pilot-induced 
oscillations  and  for  optimizing  design  alternatives  such  as  seat  loca¬ 
tion,  orientation  and  suspension  parameters.  Flight  control  system 
designers  can  make  use  of  BIODYN-80  output  in  optimizing  vehicle/  air¬ 
crew  ride  qualities  and  visual  performance  effects,  possibly  incorporat¬ 
ing  anti-vibration  devices  to  improve  the  design. 

One  further  application  of  BIODYN-80  deserves  special  mention.  Its 
vibration-input  to  blodynamlc  parameter-output  transfer  functions  are 
ideally  suited  as  input  to  PIVIB,  another  software  package  which  relates 
pilot  tracking  performance  to  the  vibration  environment  (Reference  7). 
PIVIB  accepts  biomechanical  transfer  functions  in  the  format  created  by 
BIODYN-80.  The  details  of  the  BI0DYN-80/PIVIB  interface  will  be  found 
in  a  later  section. 


The  remainder  of  this  report  details  the  use  of  BIODYN-80,  and  des¬ 
cribes  the  model  and  the  equations  comprising  it.  Detailed  instructions 
for  the  creation  of  the  required  input  files,  and  a  complete  example 
problem  are  also  included. 
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SOFTWARE  OVERVIEW 


Figure  1  is  a  functional  block  diagram  description  of  the  ele¬ 
ments  in  BIODYN  and  its  interface  with  PIVIB.  The  BIODYN-80  package  is 
composed  of  three  programs.  The  first,  called  CREATE,*  interactively 
sets  up  the  two  input  files  used  by  BIODYN.  The  second  program  called 
BIODYN,  is  the  actual  "number  cruncher,"  which  structures  and  solves  the 
biomechanical  equations  and  computes  the  desired  transfer  functions. 
The  third  program  called  PLOT,  accepts  the  file  of  BIODYN  transfer  func¬ 
tions,  prints  selected  ones  in  a  form  readily  comprehended  by  the  user, 
and  prepares  "quick  plot"  Bode  plots  on  the  line  printer,  to  facilitate 
a  visual  interpretation  of  the  transfer  function  information.  Both 
BIODYN  and  PLOT  are  designed  primarily  as  batch  programs  while,  CREATE 
permits  conversational  user  interaction  in  structuring  input  data. 

A  subsequent  link  in  this  series  of  programs  is  PIVIB.  It  Is  a 
large  batch  program  with  three  modules.  The  first,  BDMOD ,  computes  the 
response  behavior  of  the  various  biomechanical  subsystems.  The  second, 
PVMOD ,  uses  the  results  of  BDMOD  and  the  BBN  optimal  control  model  to 
estimate  pilot  tracking  performance  within  the  vibration  environment. 
The  final  module,  VEXEC,  provides  the  top  level  communication  interface 
between  BDMOD  and  PVMOD,  and  performs  no  actual  computation.  The  BDMOD 
module  expects  biodynamic  transfer  functions  in  a  format  generated  by 
BIODYN.  Only  the  relationships  of  BIODYN-80  to  PIVIB  are  described 
herein,  because  PIVIB  is  run  separately  and  has  its  own  User's  Manual 
(Reference  7). 

Figure  2  presents  a  flow  chart  covering  the  use  of  BIODYN  and  PIVIB 
in  a  given  session.  The  detailed  description  of  the  various  steps  for 
BIODYN  is  found  in  Section  D.  Note  that  the  flow  of  execution  can  be 


*The  CREATE  program  should  not  be  confused  with  INTERCOM'S  EDITOR 
subcommand  CREATE;  it  has  been  given  the  permanent  file  name  EXECRT  to 
distinguish  it  from  the  latter. 
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used  to  solve  a  single  problem,  by  submitting  a  single  batch  request,  or 
to  Iterate  on  a  design  by  submitting  a  number  of  batch  requests  using 
the  same  CHOICES  file  and  slightly  altered  PARAMETER  files  (heavy  lines 
In  Figure  2). 

BIODYN-80  uses  three  separate  files  in  performing  Its  computations. 
The  PARAMETER  file  contains  the  set  of  96  parameters  used  to  define  the 
specific  pllot/posture/di8play /vibration  characteristics.  Appendix  A 
contains  a  complete  description  of  each  parameter,  Including  definition, 
mnemonic,  nominal  value,  recommended  range  of  values  and  reference 
(where  available).  CREATE  is  used  to  assemble  this  file  and  can  modify 
an  existing  file  or  produce  an  entirely  new  file.  The  CHOICES  file 
contains  the  list  of  desired  transfer  functions  to  be  computed  and 
output  by  BIODYN,  as  well  as  directives  for  producing  the  line  printer 
Bode  plots.  Again,  this  CHOICES  file  is  assembled  by  CREATE.  Finally, 
the  TF  file  Is  used  to  store  the  resulting  transfer  functions  output  by 
BIODYN,  and  is  read  by  the  PLOT  routine  for  generating  Bode  plots. 

PIVIB  employs  a  single  large  file  to  direct  Its  flow  of  execution. 
This  file  defines  the  vibration  environment,  biomechanical  transfer 
functions,  tracking  dynamics,  tracking  performance  requirements,  and 
pilot  limitations  (bandwidths,  time  delays,  etc.).  Currently,  this  file 
Is  assembled  In  the  editor,  using  output  from  BIODYN-80  if  desired. 

BIODYN-80  and  each  of  its  predecessors  were  developed  on  the  Tym- 
share,  Inc.,  PDP-10  computers.  BIODYN-80  has  been  adapted  to  the  CDC 
"Intercom"  System  on  the  CDC  6600  or  CYBER  175  at  WPAFB  In  order  to 
increase  Its  availability  to  Air  Force  users,  and  to  Interact  with 
PIVIB,  which  is  also  operable  on  the  WPAFB  CDC  computer.  The  details  of 
this  manual  will  address  its  use  on  the  CDC  machine;  an  example  of  a 
Tymshare  session  is  given  in  Appendix  D.  Throughout  the  manual,  how¬ 
ever,  it  is  assumed  that  the  user  is  familiar  with  the  WPAFB  CYBER  175, 
and  in  particular  has  experience  with  the  INTERCOM  operating  system.  If 
not,  the  user  should  read  References  8  and  9  first. 
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MODEL  DESCRIPTION 


Three  distinct  subsystem  models  are  included  in  BI0DYN-80.  They  ere 
described  individually  below. 

A.  BIOMECHANICAL  MODEL 

Figure  3  (updated  from  Figure  2  of  Reference  3)  presents  the 
biomechanical  model  and  defines  many  of  the  necessary  parameters  that 
describe  the  nominal  (or  trim)  situation.  It  utilizes  an  "isomorphic, " 
or  body-mimicking  representation,  of  the  major  body  segments  in  their 
orientations,  simplified  to  a  minimum  number  of  lumped  parameter 
equivalents.  The  biomechanical  features  Include: 

•  Semisupine  torso;  sliding  hip  plus  rocking  chest 
supported  on  a  compliant  buttocks/seat. 

•  Head  bobbing  on  an  articulated  neck  with  passive 
compliance,  or  active  neuromuscular  system. 

•  Upper  arm  and  forearm  links  plus  grip-interface  com¬ 
pliance,  driven  by  an  active  neuromuscular  system. 

•  Arm-rest  restraints  (optional). 

•  Stick  "feel  system"  dynamics  from  zero  to  infinite 
stiffness,  and  any  angle  of  stick  or  grip. 

The  simplified  torso  model  was  derived  to  describe  the  dominant  mo¬ 
tions  of  the  head  and  arm  elements;  the  "pin  joint"  node  between  upper 
and  lower  torso  segments  is  not  meant  to  represent  any  physical  feature. 
In  practice,  the  masses  and  inertias  are  obtained  from  tabulated  bio¬ 
mechanical  and  anthropometric  data  for  the  appropriate  sized  person 
(e.g.,  Reference  11),  the  postural  angles  are  based  on  the  actual 
situation  (preferably  via  a  side-view  photo),  and  the  spring  forces  and 
damping  coefficients  are  fitted  to  data  or  taken  from  other  sources 
(e.g..  References  3  and  12). 
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a)  Basic  Elements 


Figure  3  Main  Biomechanical  F.lements 


B.  LIMB  NEUROMUSCULAR  MODEL 


The  active  neuromuscular  system  noted  on  Figure  3  Is  a  schematic 
representation  of  the  net  effect  of  complex  agonlst/antagonlat  muscle 
pairs  controlling  the  upper  an  or  head  based  on  the  work  summarized  In 
Reference  10.  An  ”NM  switch"  Is  defined  In  the  PARAMETER  file  which 
causes  this  neuromuscular  model  to  control  the  limb  (NM  -  0)  or  the  head 
(NM  ■  1).  A  linearized  representation  of  the  limb  neuromuscular  model 
Is  shown  In  Figure  4,  while  the  head  neuromuscular  model  is  depicted  In 
Figure  5.  This  model  relates  the  action  of  the  muscle  pairs  to  the 
effective  (spindle)  sensors  of  muscle  length  and  force  as  well  as 
proprioception  from  the  stick  grip  interface  (in  the  case  of  the  limb 
neuromuscular  model)  or  the  head-neck  Interface  (in  the  case  of  the  head 
neuromuscular  model),  thus  closing  the  receptor-CNS-effector  loop. 

Unless  neuromuscular  properties  are  being  Investigated,  It  Is  recom¬ 
mended  that  the  typical  values  of  the  parameters  shown  in  Table  A-l  be 
used.  These  are  representative  of  a  normal  person's  arm-hand  or  head- 
neck  system,  and  generally  yield  reasonably  damped  neuromuscular  servo 
properties.  The  neuromuscular  parameters  listed  In  Appendix  A  are 
characteristic  of  the  largest  muscles  in  the  body  (e.g.,  the  legs),  but 
experience  has  shown  that  the  dynamic  properties  (torque/lnertla  ratios, 
damping  ratios,  natural  frequencies,  etc.)  are  about  the  same  for  all 
postural  muscle  pairs.  Here,  an  empirical  scale  factor  S  (>  1.0)  is 
used  to  scale  the  normalized  muscle  to  a  particular  configuration,  as 
though  the  muscle  acted  normal  to  the  upper  arm  c.g.  A  more  detailed 
description  of  the  neuromuscular  model  is  found  in  Reference  10. 

C.  HEGK/HEAD/EYE  MODELS 

Figure  6  illustrates  the  basic  elements  involved  in  the  neck/head/ 
eye  model.  Details  involving  validation  and  example  use  of  this  subsys¬ 
tem  model  are  found  in  Reference  6.  The  task  is  to  keep  the  eye  fixated 
on  the  target,  i.e.,  null  the  relative  eye  (point  of  regard)  deviation 
(RED)  at  the  display.  The  moving  base  can  produce  image  error  distur¬ 
bances,  both  from  induced  head  rotation  and  translation  as  well  as 
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Simplified  Head  Pointing  Neuromuscular  System 
(When  Option  NM  »  1) 


s)(2  -  rus)  |  Angular  Image  Motion  Feedback 


I 

i 


I 


Figure  6  Block  Diagram  for  Head,  F.ve,  Vestibular  Models  and  Display  Kinematics 


target  motions.  The  "vestibulo-ocular  reflex”  is  a  "crossfeed"  (or 
"feedforward”)  that  rotates  the  eye  oppositely  to  the  head  to  compensate 
for  head  rotation,  i.e.,  to  maintain  approximate  inertial  eye  fixation. 

The  "Fixation  Reflex"  or  feedback  tracking  loop  involves  the  sub¬ 
ject's  efforts  to  null  the  image  error  0j  by  compensatory  eye  movements. 
The  "Target  Pursuit"  path  (shown  dashed)  models  the  subject's  "feedfor¬ 
ward"  operations  on  the  perceived  absolute  target  motions  in  inertial 
space  (as  distinct  from  the  image  error  motions).  For  highly  predict¬ 
able  and  perceivable  target  motions  the  Target  Pursuit  path  is  capable 
of  greatly  reducing  the  image  errors,  but  it  is  not  involved  in  most 
vibration  cases  because  of  insufficient  bandwidth. 

The  linearized  model  formulation  assumes  that  the  target  remains 
within  the  foveal  area  (3-4  deg  field),  and  that  angular  velocities  are 
small  (i.e.,  less  than  20  deg/sec.)  such  that  saccades  are  rare  and  can 
be  ignored.  Additionally,  the  vibration-induced  motions  will  always  be 
small  enough  to  permit  linearization  of  all  angular  functions  and  allow 
a  quasi-linear  representation  of  the  dynamic  elements  about  each  "oper¬ 
ating  point"  (posture,  view  geometry,  frequency).  The  operating  point 
angles  can  have  any  values  within  human  limits  for  a  seated  pilot  or 
crewman. 

The  three  subsystem  models  described  above  contain  parameters  which 
must  be  assigned  values  via  the  PARAMETER  file,  in  order  to  describe  the 
desired  problem  to  BIODYN.  These  parameters  are  defined  in  Appendix  A. 
The  available  input  (forcing  function)  variables  and  output  (response) 
variables  are  given  in  Appendix  B.  Two  different  input  deflection  vari¬ 
ables  are  available  for  the  x  and  z  axes.  BZP  and  BXP  should  be  used 
when  absolute  ouputs  (referenced  to  inertial  spacr  '  are  desired;  DZP  and 
DXP  are  used  when  outputs  relative  to  the  moving  platform  are  desired. 
The  equations  of  motion  are  in  Appendix  C.  Background  on  their  deriva¬ 
tion  and  development  are  given  in  a  series  of  prior  reports  (References 
1-3,  6,  10). 
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USING  BIODYS— 80 


Three  job  steps  an  required  for  a  complete  run  of  BIODYN-80.  The 
CREATE  program  is  used  interactively  via  a  remote  terminal,  to  assemble 
the  PARAMETER  and  CHOICES  files.  These  files  are  then  read  by  BIODYN  and 
the  transfer  functions  computed,  stored  on  the  TF  file,  and  listed  at  a 
line  printer.  The  TF  file  is  in  turn  read  by  the  PLOT  program  to  pro¬ 
duce  any  Bode  plots  desired.  Each  of  these  steps  is  described  in 
detail.  Note  that  all  user  Interaction  with  the  computer  assumes  the 
user  is  logged  in  to  CDC's  INTERCOM.  It  is  assumed  that  the  user  has 
access  to  Intercom  set  "AFMLH  where  the  program  and  its  files  reside. 
For  other  setname  groups,  please  consult  your  Computer  Center  Represen¬ 
tative  to  get  on  the  AFML  list,  or  for  instructions  on  how  to  use  a 
"modified  Job  stream  listing"  (Reference  8). 

A.  CREATE 

Because  of  the  large  number  of  parameters  required  by  BIODYN-80,  it 
was  decided  that  a  dedicated  Interactive  "editor"  for  creating  parameter 
files  was  needed  to  produce  rational,  error-free  input  files,  and  there¬ 
by  avoid  aborting  during  the  subsequent,  and  expensive,  batch  run  of 
BIOOYN  (for  example,  due  to  a  missing  decimal  point).  Furthermore,  this 
approach  allows  the  maintenance  of  a  "catalog”  of  PARAMETER  files  for 
various  postures  and  situations,  from  which  the  most  appropriate  one  can 
be  easily  selected  and  modified  for  a  particular  situation.  The  CREATE 
program  serves  in  this  regard  by  allowing  a  user  to  interactively  assem¬ 
ble  new  parameter  files  or  modify  existing  ones.  CREATE  also  performs 
various  data  checks  and  identifies  common  parameter  entry  errors.  This 
section  details  the  use  of  the  CREATE  program. 

1 .  Entering  CREATE 

CREATE  makes  use  of  five  files,  as  defined  by  the  PROGRAM  card: 

■  INPUT  file  (the  user's  terminal) 
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TAPE4 


TAPE  5 


OUTPUT  file  (the  user's  terminal) 


T 


TAPE7  »  CHOICES  file  (if  new  one  is  created) 

TAPE8  -  PARAMETER  file  (new  or  modified) 

TAPE 20  -  PARAMETER  file  (existing) 

The  PROGRAM  card  automatically  assigns  the  INPUT  and  OUTPUT  files  to  the 
appropriate  tape  units.  If  an  existing  PARAMETER  file  is  to  be  ac¬ 
cessed,  it  must  be  assigned  prior  to  executing  CREATE.  The  job  stream 
to  be  input  by  the  user  looks  like  this  (detailed  examples  are  given 
later,  in  Appendices  D  and  E): 

REQUEST, TAPE 7, *PF 
REQUEST, TAPE8,*PF 
ATTACH, TAPE20,pfn^ 

ATTACH, EXECRT 
EXECRT 


(pfn^  is  the  PARAMETER  file  name 
selected  by  the  user  from  the 
permanent  file  catalog  within 
the  set name  AFML  group) 

I 


Now  the  user  is  at  the  entry  level  of  CREATE.  The  first  issue 
addressed  is  creating  the  PARAMETER  file,  or  modifying  an  existing 
PARAMETER  file.  The  computer's  first  prompt 

NEW  FILE? 

asks  the  user  to  identify  which  option  is  desired.  Answers  are  YES  or 
NO;  any  other  response  will  cause  the  query  to  be  repeated. 

2.  Modifying  a  PARAMETER  File 

If  the  answer  is  NO,  the  values  stored  on  the  file  assigned  to 
TAPE20  are  read.  The  user  may  list  those  values  by  responding  to  the 
next  prompt 

LISTING  DESIRED? 
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with  YES.  A  NO  response  will  skip  over  the  listing  option.  Any  other 
response  will  csuse  repeat  of  the  question.  The  listing  produced  can 
either  be  long  or  short.  The  long  listing  provides  parameter  defini¬ 
tion,  mnemonic,  units  and  current  value,  while  the  short  listing  prints 
only  values,  grouped  according  to  each  element  (see  Table  A-2,  Appen¬ 
dix  A).  Examples  of  each  are  given  in  Figures  7  and  8,  respectively. 

Response  to  the  prompt: 

(LO)NG  OR  (SH)ORT  LISTING? 

determines  which  listing  to  generate.  A  response  other  than  LO  or  SH 
will  cause  repeat  of  the  question. 

Values  are  now  ready  for  change,  as  Indicated  by  the  message: 

INPUT  MNEMONIC  (5  CHARS)  AND  VALUE  FOR  EACH  CHANGE, 

TO  EXIT,  TYPE  XXX. 


To  change  any  desired  parameter  value,  its  mnemonic  from  the  second 
column  of  the  definitions  in  Appendix  A  is  input.  Mnemonic  must  be  left 
Justified  within  a  5  character  field  followed  by  a  comma,  followed  by 
the  new  value.  The  mnemonic  is  then  checked  against  the  internal  list 
of  permissible  mnemonics;  e.g.,  if  ABDCE  is  not  found,  the  message: 

MNEMONIC  ABCDE  NOT  PERMISSIBLE.  PLEASE  REINPUT. 

appears  and  the  user  should  try  again.  Likely  errors  include  typos,  and 
forgetting  to  pad  the  5  character  field  with  trailing  blanks. 

Once  a  correct  mnemonic  is  identified,  its  new  value  is  compared 
with  the  recommended  range  of  values  for  that  parameter.  If  the  value 
is  not  within  that  range,  a  warning  message  appears: 

WARNING  -  RECOMMENDED  RANGE  FOR  THIS  PARAMETER  IS 
0.5  TO  3.5  GS 

DO  YOU  WANT  TO  CHANGE  THE  VALUE? 
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STD  CREWMAN »  HANDS  ON  KNEES 


PARAMETER  Ml:  i  INI 

T  1  ON 

MNEMONIC 

UNITS 

VALUE 

LOWER  BODY  MASS 

XMB 

KG 

.1.4,0000 

i.  •  B  0  D  Y  3  1C  r'i  T  i  U  S I  I 1 0  N 

DAMPER 

BB 

N/M/S 

315 • 000 

L  .  .BODY  SEAT  CUSHION 

ST  ILLNESS 

XKB 

N/M 

49360.0 

SERIES  SPRING  GRAD, - 

LOWER  BODY 

GK 1 

N/M 

119320. 

SEAT  ill.  I  ANGLE 

THE 

DEG 

13.0000 

fE:i  M  1  'E  E  i  Ji.iAi..  ACCEL  . 

G 

M/S2 

9.80000 

TORSO  MACS 

XMT 

KG 

1.8,0000 

TORSO  INERTIA 

zrr 

KGM2 

.800000 

TORSO  DAMPING 

BTB 

NM/R/S 

40.0000 

TORSO  S  i  !  I:  ENESS 

XKTB 

NM/R 

500 . 000 

ANGLE'  OE  X i ...  T  FROM  VE 

RT  LCAI 

THT 

DEG 

-15.0000 

TORSO  CG  TO  HIPS  PIVOT  LENGTH 

XLT 

h 

. 150000 

ANGI  E  OF  XI  TN 

THIN 

DEG 

-5.00000 

TORSO  GO  TO  NECK/ TORSO  PIVOT 

XLTN 

M 

.300000 

ANG!  \"  TK  XL S 

THS 

DEG 

-5.00000 

TORSO  T.G  TO  SHOULDER 

LENGTH 

XLS 

M 

.300000 

NEC  K  M  A  S  S 

XMN 

KG 

1,24000 

NECK  INERTIA 

ZIN 

KGM2 

, 240000E -02 

NECK/ TORSO  DAMPER 

BNT 

NM/R/S 

. 100000 

NECK /TORSO  STIFFNESS 

XKNT 

NM/R 

50.0000 

ANGLE  OF  XI  N 

THN 

DEG 

-30.0000 

NECK  LENGTH 

XLN 

h 

. 124440 

NECK  r B  TO  NECK.  TORS 

0  PIVOT 

XLN1 

M 

. 622200E-01 

HEAD  MASS 

XMH 

KG 

3.10000 

HEAD  INF  R’I  V  A 

ZIH 

KGM2 

,  303000E -0.1 

HEAD/NECK  DAMPER 

BHN 

Nh/R/S 

.126000 

HE AI) /NECK  STIFFNESS 

XKHN 

NM/R 

15.0000 

HE AD /NECK  COMPLIANCE 

CH 

- 

1.00000 

ANGLE  OF  /LEI 

TH 

DEG 

60.0000 

HEAD  CG  iO  Hi  All/ NIH'K 

PIVOT  I..EN 

XLH 

M 

. 250000E-01 

LINE  OE  SIGHT  ANGLE 

THV 

DEG 

-30.0000 

VIEWING  DISTANCE 

VD 

M 

.685300 

UPPER  ARM  MASS 

XM1 

KG 

1 . 37200 

UPPER  ARM  INERTIA 

Z1 

KGM2 

. 120000E-01 

UPPER  ARM  ANGLE 

ri 

DEG 

15.0000 

UPPER  ARM  LENGTH 

D1 

M 

. 290000 

UPPER  ARM  CG 

Ri 

M/M  • 

.440000 

LOWER  ARM  MASS 

XM2 

KG 

1.01700 

LOWER  ARM  INERTIA 

7.2 

KGM2 

, 152000E-01 

ELBOW  ANGLE 

TE 

DEG 

115.000 

LOWER  ARM  LENGTH 

D2 

M 

.305000 

LOWER  ARM  CG 

R2. 

M/M 

.500000 

GR IP  INTERLACE  ANGLE 

TIJ 

DEG 

0, 

GRIP  INTERFACE  TIME 

CONSTANT 

TIF 

S 

.100000 

GRIP  INTERFACE  COMF1.. 

I A  MCE 

Cl 

M/N 

(continued) 

1.00000 

Figure  7  Example  Long  Format  File  Listing 
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■  ’  i  :  l'1.  V-  • 

XMS 

EG 

4 .  OOOOO 

■Ml  ■ 

BS 

N/M  -  S 

? 4 0.000 

:  |  !  !■  i.  ;|-  .  1  ( !  1  i  H  ! 

XKS 

N/M 

1 0000 .0 

•  1  ;  1  ;  :  :  It!  i  !  ’.tJ!  :  ; 

nKiVilr  1  l::.l  v 

CKS 

.... 

I  . OOOOO 

rfi  iiU-ii  !  !'.i : :  , i  s-*f.»  i  i  i. 

i i v  ro  nx 

SX 

N/G 

0  ♦ 

Hi  :  HWl::  hill  S|-  r>  •  >  i  i  1  v 

1.  1  f  TO  NZ 

sz 

N  /  b 

0. 

.  „  '1  .|  i: 

THE 

DEG 

?0 *  0000 

■■■  •  i  -  »  i  !-  rn-.i  f  !  ■! 

xlc: 

M 

. 55 00 00 

MUlRiif  Si  Mi  1 

i  i  Hi  i  uk 

XKSC 

/M 

*  1.0 0 000 IF  ~  0 

t1-.;  ',  t-  i '  S  't  iril  V  K  »  i- 

ORMAE > 

BAR 

N/M/S 

0 . 

'V-  f  '  t  i  1.  !•  {■  Nf  /• 

< NOR HAL ) 

XKAR 

N/M 

0 ,. 

'll '  E ,  .  in  l hr  !'  P  ■’  T 

A  NO  > 

BA  f 

N/M/S 

0  v 

■  v  n  ;  T  i  i-i  N!  '••• 

i  !  A  N  Ci ) 

XKAT 

N .  M 

0  ♦ 

'  i-.  !  Li  !  1  i  }  •  H  1“  H  f 

ti  !  :  i  (  AiA  ~ 

XLE'R 

M 

0  * 

■  K  :'V  . ') }•■.’ M  u)!-  |  •..!'•(  )’ 

i  :  1  '  i*l  1  :  i;:  0  1 

ARMS' 

.... 

o  „ 

•  U  !  :  T  }  i  !•  OK  A ! 

ANT:  1--  A  *  Tl  * 

NM 

- 

()  v 

jVRPmI  .  '.!•  !  !':!•,  H  S 

;  1  1  1  r:  HE  NTS 

s  :i: 

GO . 0000 

I  0  !  ?V  1  ,  1  1 

roc 

S 

9  ,  y  j{.  .  ,  j..  . 

■!  :'■*  :  i ;  •  1  '  ’  .  "i  T  |  :  if 

Oi  l  iT  r 

res 

.  S 9 00 00 IF  -0 

Pill!  .  '!  .  .  !■ 

XT-  A  A 

N/N 

1  » OOOOO 

:  •  i  i  i  -l  ,i  !  ;  .1  |-  .  ,  i  ft  v 

ZA 

.... 

.  3 0 0 0  0 0 

■1.1  i  i.ir  f,l  f  I-'!.  )!  !l  Ml!  ;■ 

r  N  T  f:i .  0  , 

WA 

R/S 

:l.  A  ,  0000 

i  i  ll  ill  -  i  i  i  1  !  i'J  1  .  A  . 

0 .. 

TAA 

s 

0 . 

M  !  :  i  |  iiUI  ■  'm?  tl  !  I ' 

BM 

N/N/S 

1  .  OOOOO 

i!  !■■  !  f.'l.i  !  r-i  M  h  Vr  r 

EM 

X  K  M 

N/M 

2 . OOOOO 

Si:  P  If.-.  !'  i  A--.  (  |.i  •;  (.  F 

HEN T  TAMPER 

BE 

N/M/S 

2 . 43000 

'■■:  !  !■■  >•.  '  i  A  ;  •  1"  1  I  1 

HI  N  ?  OR  Ail 

XKE 

N/M 

40 . 0000 

!"■■  ,  1;  J  P.-thpi  !'■ 

BT 

N/M/S 

0.. 

i  r  'll  •  Ml  .  H  ' 

XKT 

N  /M 

80 . 0000 

i  "III.  :  rsi/iil"  i  A  A  i'N 

XI\SP 

N  /  h 

5 .  OOOOO 

i  •  - ! : !  :  !  r:  Ail  r  1  f  i ! 

i  i)  NS  TAN  I 

ISP 

s 

» 909 0001:. -0 

A  •  1 .1 ! • '  !  !  .  .!■■  i  !  hi- 

0 NS T ANT 

TPS 

C; 

0, 

•••'  i  ■•'•'!  !  'll  !  T  !'  1  1: 

1 1  IT  I  .  C  . 

TSS 

•“) 

0, 

; 1  ■  i  !.!  :■  :  i  .ir  i  M 

XKG 

N/N 

. 500000 

1  "  1  '  •  :P  i  i  III-  At; 

r!  S  i  (")  M  I 

TG 

i:> 

*  550000E.  ~"0 

I':':  -  ■  !  i  !  i  :  '  I  :-  1  !  :  1  1 

i IMP  OONST 

TZ 

s 

\J  * 

■  ;  1  '  Ai  •  Mr 

s  f  AN  | 

TP 

c;; 

0, 

•'  •  i  '  »  :.'l  i 

;  .1 

XKBE 

R  /  R 

7 . 07950 

XKRE 

R/R  /S 

0. 

i  ION  I..OOP 

TO 

S 

>  450000F- 0 

■  '  ’  .  :  1  '  •  1  1 

'•ill  EE  El. • 

AI.  PH  A 

R/S 

. 300000 

.  ! 

■  i 

XKP 

R/R 

t 

■  :  H  SI  NSR 

XKIE 

.... 

.  670000 

!  I E  SEN SR 

X  l\  v !: 

- 

, 100000 

re? 

s 

. 1 000 GOT  0 

■  ■  . 

■ '  ; 

TEM 

s 

4 ] OOOOO 

:  AM  i 

n.H 

s 

4 1 2 50 OOF -0 

. iif*  s  i'i'*  . 

l  HOP! 

OE 

J  4  0  0  0  0  0 

i  : 

ZE 

.... 

. 650000 

>  .  t  '  •  . 

!  ■  : 

WP 

R/S 

3  i.  A  ..  230 

Figure  7 

(Concluded) 
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PARAMETERS  FOR 


SEATED  CREWMAN  -  VIEWING  DISTANCE  •-  9  CM 


LOWER  BODY 

14.000 

3 1 5 . 00 

49360. 

.  1 1932.E+06 

13.000 

9.8000 

TORSO 

18 . 000 

.80000 

40 . 000 

500 . 00 

- 15.000 

.15000 

-5 . 0000 

. 30000 

-5 . 000 0 

. 30000 

NECK 

1  ,  .'400 

. 24000F  02 

.  10000 

SO . 000 

-JO .000 

. 12444 

. 62220E-01 

HEAD/DI SPLAY 

VIEWING 

3 . 1 000 

, 30300E-0 1 

. 12600 

1 5 . 000 

1 . 0000 

AO . 000 

.25000E-01 

-30.000 

,  9000 OE -0.1. 

ARM  <  UPPER t 

LOWER) 

1 . 3720 

.  12000E-01 

15.000 

. 29000 

.44000 

1 .0170 

, 15200E-01 

115.00 

, 30500 

,  5  0  0 0 0 

GRIP  INTERFACE 

0. 

. 1 0000 

1.0000 

STICK 

4.0000 

240.00 

10000. 

1 . 0000 

0. 

0 , 

90.000 

. 55000 

. 10000E-01 

ARM  REST 

0. 

0 . 

0 . 

0. 

0 . 

o  • 

NEU  ROMUSCIJLA 

R  SYSTEM 

0  50.000 

.  20000E -01 

. 90900E-01 

. S9000E-01 

1 , 0000 

.80000 

16.000 

0. 

1.0000 

2 . 0000 

2.4300 

40 . 000 

0 . 

8 0 . 0 0 0 

5 . 0000 

.  90900E -0.1 

0. 

0. 

.50000 

.55000E-01 

0. 

0. 

IMAGE  F I X AT 1 ON/ VEST I BULO 

-OCULAR  SERVO 

7 . 0795 

0 . 

. 45000E-0 1 

. 30000 

0 . 

.67000 

.  10000 

. 10000E-01 

. 10000 

,  12500E-01 

1.0000 

. 65000 

3 16,23 

Figure  8  Example  Short  Format  PARAMETER  File  Listing 
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If  a  mistake  has  been  made  In  typing  the  new  value,  a  YES  answer  will 
allow  retyping  of  both  the  mnemonic  and  Its  corrected  value.  If,  how¬ 
ever,  the  value  Is  correct,  even  If  out  of  range,  e  NO  answer  will  cauae 
the  value  to  be  stored  as  it  stands.  Any  response  other  then  YES  or  NO 
will  cause  the  warning  aassega  to  be  repeated. 

The  range  of  parameters  given  In  Appendix  A  la  as  wide  as  can  be 
allowed  for  typical  applications  of  BI0DYN-80  (e.g.,  arm  joint  limits). 
Values  outside  this  range  are  used  at  the  user's  risk,  and  precise  pos¬ 
tural  angles  for  the  given  case  should  be  carefully  checked  at  this 
point.  (It  la  good  practice  to  check  your  lengths  and  angles  by  drawing 
a  stick  figure  to  scale,  using  the  desired  values.) 

When  all  the  desired  changes  have  been  made,  an  Input  of  XXX  causes 
control  to  leave  the  modify  mode.  An  opportunity  to  change  the  case 
title  is  then  presented  as 

TITLE  IS  STD  CREWMAN,  HANDS  ON  KNEES 

CHANGE  DESIRED? 

A  YES  answer  causes  the  message 

NEW  TITLE: 


to  appear  and  a  new  title  (60  characters  maximum)  can  be  typed.  A  NO 
response  skips  this  section;  other  responses  cause  the  prompt  to  repeat. 

After  all  desired  changes  have  been  made,  the  user  is  once  again 
given  the  opportunity  to  list  the  file,  using  either  the  LOng  or  SHort 
format,  as  described  above.  At  this  point  the  computer  writes  the  new 
values  to  the  PARAMETER  output  file,  TAPE8. 

3.  Assembling  a  New  PARAMETER  File 

If  most  of  the  parameters  are  new,  then  instead  of  modifying  an 
existing  file  the  initial  prompt 
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NEW  FILE? 


Is  answered  YES.  Every  parameter  value  must  then  be  input.  The  message 

INPUT  VALUE  FOR  EACH  MNEMONIC.  FOR  FURTHER  EXPLANATION,  TYPE? 

Initiates  this  process.  For  each  parameter,  the  mnemonic  is  printed  as, 
for  example,  Viewing  Distance,  VD: 


VD  - 

and  the  value  is  accepted  in  floating  point  format.  If  the  user  is 
uncertain  as  to  the  range  of  appropriate  values,  an  expanded  prompt  may 
be  requested  by  typing  ”?”.  The  expanded  prompt  for  VD  is,  for  example: 

VD  -  T_ 

VD  VIEWING  DISTANCE 
RECOMMENDED  RANGE  -  0.5  TO  3.0  M 
VD  - 

the  user  can  then  enter  the  desired  value. 

As  each  value  is  read,  it  is  checked  against  the  recommended  range 
of  values,  in  order  to  weed  out  input  errors  due  to  wrong  sign  or  incor¬ 
rect  decimal  placement.  If  the  value  is  not  within  range,  the  warning 
message  given  previously  appears.  The  user  always  has  the  option  of 
modifying  or  saving  any  value  typed. 

This  process  continues  until  all  ?6  parameters  have  been  entered. 
The  case  title  is  then  added  as  a  response  to 

INPUT  TITLE  FOR  THIS  CASE: 

At  this  point,  a  LOng  or  SHort  format  listing  may  be  generated  as  dis¬ 
cussed  above.  The  file  is  then  written  by  the  computer  to  the  PARAMETER 
file  designated  TAPE8. 


4.  Assembling  the  CHOICES  File 

The  next  section  of  CREATE  is  used  to  assemble  a  new  CHOICES  file 
which  specifies  the  transfer  functions  to  be  output.  Note  that  there  is 
no  provision  for  modifying  an  existing  CHOICES  file;  any  changes  to  an 
existing  CHOICES  file  must  be  made  by  creating  an  entirely  new  file. 
The  prompt 

NEW  CHOICES  FILE? 


has  two  responses:  a  NO  allows  the  user  to  skip  this  entire  section  if 
an  existing  CHOICES  file  will  be  used  as  input  to  BIODYN;  a  YES  response 
will  cause  prompting  for  the  CHOICES  file  components  as  described  below. 
Any  other  response  will  cause  the  prompt  to  be  repeated. 

A  YES  response  to  the  next  prompt 
BIODYN  TFS  DESIRED  FOR  PIVIB? 


will  automatically  generate  the  ten  transfer  functions  to  be  Included  in 
the  PIVIB  input  file.  This  is  mandatory  if  Interaction  with  PIVIB  is 
anticipated.  The  response  and  forcing  function  variables  involved  are 
shown  later,  in  Figure  12.  A  NO  answer  skips  this  section;  other 
answers  cause  repeat  of  the  question. 

All  other  desired  transfer  functions  are  entered  by  the  user  in  re¬ 
sponse  to  the  following  computer-specified  format: 


TRANSFER  FUNCTION  INPUT: 

FIRST  LINE  -  RESPONSE  MNEMONIC,  FORCING  FUNCTION  MNEMONIC 
(AAA,  AAA);  ENTER  XXX  TO  STOP 
SECOND  LINE  -  PLOTTING  INFORMATION,  5  ITEMS: 

BODE  LOWER  FREQ.  LIMIT 

BODE  UPPER  FREQ.  LIMIT 

BODE  UPPER  PHASE  LIMIT  (0.  DEFAULTS  TO  200.) 

BODE  LOWER  PHASE  LIMIT  (0.  DEFAULTS  TO  -400.) 

LIST  (1.  TO  LIST  TABLE,  0.  FOR  NO  LIST) 

IF  NO  PLOT  DESIRED,  ENTER  0.  FOR  ALL  ITEMS 
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Note  that  the  parameters  for  each  transfer  function  are  entered  on  two 
consecutive  lines,  each  Item  separated  by  a  comma,  each  line  followed  by 
CR. 


The  mnemonics  requested  are  three-character  codes  as  defined  In 
Appendix  B.  The  program  checks  each  mnemonic  against  the  list  of 
acceptable  codes.  If  a  match  is  not  made,  the  error  message 

XYZ  NOT  PERMISSIBLE,  PLEASE  REINPUT, 
is  given  and  the  entire  first  line  must  be  retyped. 

If  Bode  plots  are  desired,  the  frequency  and  phase  limits  must  be 
entered;  if  not,  zeros  are  assumed.  The  frequency  limits  must  bound  a 
frequency  range  of  no  more  than  3  decades.  If  frequencies  with  an 
unacceptable  range  are  entered,  the  message 

MAX  FREQUENCY  RANGE  IS  3  DECADES 

PLEASE  REINPUT  ENTIRE  LINE 


appears  and,  all  five  items  for  that  transfer  function  must  be  re-typed. 
The  phase  limits  are  unbounded,  but  entering  zeros  will  cause  default  to 
+200.0  deg  and  -400.0  deg  for  the  upper  and  lower  limits,  respectively. 
If  a  tabular  listing,  including  frequency,  amplitude,  and  phase  at  20 
increments /decade  intervals  is  desired,  a  1.0  should  be  entered  as  the 
fifth  item;  otherwise  enter  0.0. 

The  CREATE  program  stops  when  an  XXX  is  read  in  the  CHOICES  file 
section.  The  two  files  created  must  be  saved  as  permanent  files  by 
giving  the  following  INTERCOM  commands: 

CATALOG ,TAPE7,pfn2,RP-10 
CATALOG ,TAPE8 , pfn3 , RP-10 

where  pfn2  is  a  new  (seven  character)  permanent  file  name  assigned  to 
the  CHOICES  file  and  pfn3  is  a  new  (seven  character)  permanent  file  name 
assigned  to  the  PARAMETER  file.  The  user  is  now  ready  to  exercise 
BI0DYN. 

The  user  can  interactively  create  a  number  of  PARAMETER  and  CHOICES 
files,  each  uniquely  named  for  subsequent  use  in  the  BI0DYN  program.  In 
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most  cases  one  same  set  of  transfer  function  choices  will  be  adequate 
for  a  variety  of  different  parameter  sets. 

B.  Bloom 

Because  of  memory  size  limitations  presently  Imposed  by  the  INTERCOM 
operating  system,  BIODYN  must  be  run  In  the  hatch  mode.  A  batch  Job  can 
be  submitted  through  INTERCOM  using  the  following  command  sequence 


EDITOR 
CREATE, S 

ABC, CM1 50000, STCSA. 

ATTACH, TAPE 7, pfn£.  User  supplies  desired 

ATTACH, TAPE8,pfn3.  file  name  for  each  pfn 

REQUEST, TAPE 19, *FF. 

ATTACH, EXEBIO. 

EXEBIO(TAPE4 .OUTPUT ,TAPE7 ,TAPE8 .TAPE19 .TAPE21 ) . 
CATALOG , TAPE 19 , pf n4 , RP-999 . 

*E0R 

SAVE.GOBIO.NOSEQ 

END 

COMMAND  -  BATCH, GOBIO, INPUT, HERE 


The  first  section  of  this  batch  request  assigns  the  PARAMETER  and 
CHOICES  files  produced  by  CREATE  to  the  tape  units  which  BIODYN  reads. 
The  BIODYN  program  Is  then  retrieved  and  executed- 

BIODYN  first  reads  the  PARAMETER  file  and  uses  It  to  set  up  the 
equations  In  the  form 


Ax  -  Bu 


where  x  is  a  vector  of  response  variables  (currently  dimensioned  48)  and 
u  is  a  vector  of  forcing  function  variables  (dimensioned  9).  The  de- 
sired  transfer  functions  (as  listed  In  the  CHOICES  file)  are  then  evalu¬ 
ated,  using  Cramer's  rule  and  several  advanced  factorization  algorithms 
In  order  to  obtain  the  first-  and  second-order  poles  and  zeros.  These 
are  all  written  to  a  file  called  TAPE 19  (the  TF  file),  as  well  as  to  the 
line  printer  (via  the  OUTPUT  file).  The  format  used  by  the  program  for 
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printing  the  transfer  functions  is  shown  in  Figure  9,  while  the  foraat 
for  storing  these  transfer  functions  on  TAPE19  is  given  in  Figure  10. 

The  TF  file  (TAPE19)  is  used  by  the  PLOT  program  for  generating  the 
quick-look  Bode  plots  on  the  line  printer.  Thus,  the  last  lines  sup¬ 
plied  by  the  user  in  the  batch  job  stream  "rewind”  this  file  and  assign 
it  a  permanent  file  name  pfn4. 

c.  not 

The  Bode  plots,  if  desired,  may  be  generated  at  the  user's  terminal 
or  written  to  a  local  file  for  later  examination,  as  the  user  wishes. 
The  user's  job  stream  looks  like  the  following: 

ATTACH .TAPE 19, pfn4.  pfn4  is  the  file  name  of 

ATTACH, EXEPLT.  the  plot  instructions 

EXEPLT ( OUTPUT , TAPE 1 9 ) 

The  output  from  this  routine  consists  of  four  items.  First,  the 
input  file  mnemonics  are  printed,  so  the  user  can  Identify  which  plots 
are  forthcoming.  Then,  for  each  transfer  function,  three  items  are 
printed.  The  transfer  function  dipoles  (first  and  second  order)  which 
have  been  cancelled  are  first,  followed  by  the  numerator  and  denominator 
of  the  transfer  function.  Then,  the  Bode  plot  is  "drawn."  Finally,  if 
a  listing  was  requested  in  the  CHOICES  file,  it  follows  the  Bode  plot. 
An  annotated  example  is  given  in  Figures  11a,  b,  and  c. 

D.  UrmPAOE  WITH  PITH 

Ten  of  the  possible  BIODYN-80  transfer  functions  can  be  used  as 
input  to  PIVIB.  These  are  listed  in  Figure  12.  The  input  file 
structure  for  PIVIB  is  quite  lengthy  and  complex,  involving  over  200 
parameters.  An  explanation  of  it  is  beyond  the  scope  of  this  user's 
manual,  and  the  interested  reader  is  referred  to  Reference  7,  the  PIVIB 
User's  Manual.  A  brief  outline  of  the  steps  used  to  generate  the  PIVIB 
input  file  is  given  below. 
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Each  transfer  function  is  to  b»*  interpreted  as 


^output 

Transfer  Function  =  */lEut -  (a) 
A 


•H*  »£)}•{■»(  »?*  )} 

•{"(■’*  Vp* 


where  l  denotes  a  product  of  first-  or  second-order  roots,  and  k  Is  the  so-called  "root  locus  gain"  (of  the  high- 
frequency  asymptote): 


High  Frequency  Gain  of  Numerator 
High  Frequency  Gain  of  Denominator 


ar.d  the  various  first-  and  second-order  poles  and  zeros  are  indicated  below.  The  transfer  function  can  also  be 
interpreted  in  Bode  format,  as: 


Transfer  Function 


{n(V.  i)}.|n(^*2£.  i)j 


where  KqL  is  the  gain  of  the  low-frequency  asymptote  ("Bode  gain"): 


Low  Frequency  Gain  of  Numerator 
Low  Frequency  Gain  of  Denominator 


1 1 -MAR-80  17106 


Date,  Time 


CASE : HEAD  ?  NECK  »  TORSO  ONLY 


*  Case  Identifier 


denominator: 


.  I0464E-12 


( 

46.334  ) 

<  100.00 

)  (  330.63 

>  - 

<  •: 

.40231 

8.3929 

.  3.7123  . 

7.5273 

)  ) 

<  ( 

.42607 

10.483 

*  4 . 4663  * 

9.4335 

)  ) 

( ( 

.13334 

17.112 

*  2.1 996 

16.970 

)  ) 

( ( 

. 97433F-01 

32.746 

*  3.1912  » 

32.590 

)  ) 

( ( 

.  11386 

44.330 

»  5.0704  » 

44.241 

)  ) 

<  < 

.63021 

316.33 

f  206.02  f 

240.70 

>  ) 

.  I  2439E4- 1  2 


NUMERATOR:  KM0/P7F' 


( 

.00000 

534.5? 

)  «.  . 00000 
) 

)  <  46.384 

)  < 

100.' 

<  ( 

.44231 

*  8.3929 

*  3.7123  * 

7.5273 

>  ) 

<  ( 

.67337 

.  17.729 

f  11.94?  * 

13.104 

)  ) 

(  ( 

-.33674 

.  35.403 

• -12.631  , 

33.075 

)  ) 

<  < 

.29926 

.  44.192 

»  1 3 . 223  f 

42,166 

)  ) 

( ( 

.63024 

.  316.83 

f  206.02  r 

240.70 

)  ) 

- .  1 1234EM  2 


V 


TRANSFER  FUNCTION  COMPONENTS: 

- High-frequency  gain 

(denominator)  Kp 

- First-order  poles  ( 1  /T) 


Second-order  poles 

((?,  “Vi  -  r2)) 

Re  Im 


Low-frequency  gain  (denominator) 

Numerator •  output /input 
<High-freq.  galn>  numerator  alone 

<<High-freq.  galn>>  for  transfer  function 
K  -  K«/Kp 

First-order  zeros  r 


■  Second-order  zeros 

«Low-freq.  gain,  transfer  fUnctlon>> 
.  <Low-freq.  gain,  numerator  alone > 


Figure  9 


Format  for  BIODYN-80  Transfer 


Function  Printouts 
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-  -.'I'lWroWlim'  Wmi  11.1*11  ■  n 


Format 


Comment 


Case  title  (15A4) 
D$$ 


NP1,  NP2,  Kp,  0.,  0., 


1/T 

1/T 


Pi 

P2 


o.. 


o.. 


0. 


(Specifies  denominator) 

(214,  G14.6,  5F7.1) 

First-order  poles  (5X,  G14.6) 


CP1*  uPl  Second-order  poles  (Cl  >  “>i) 

CP2,  o>p2  (5X,  2F14.6) 


xxx/xxx 

NZ1 ,  NZ2,  Kz,  WI,  MF,  XU,  XL,  XP 

1/TZ1 

i/tz2 


4*1 1  '°*1 
4*2*  “*2 


XXX/XXX 


First  transfer  function 
numerator 

(214,  G14.6,  5F7.1) 
First-order  zeros  (5X,  G14.6) 


Second-order  zeros  (c2»  “2) 
(5X,  2G14.6) 

Second  transfer  function 
numerator 


where 


NPl  “  No.  first-order  poles 
NP2  ■  No.  second-order  poles 
Kp  -  High-frequency  denominator 
gain 

NZ1  •  No.  first-order  zeros 
NZ2  ■  No.  second-order  zeros 


K  _  High-frequency  numerator  gain 
WI  ■  Bode  plot  lower  frequency  limit 

WF  ■  Bode  plot  upper  frequency  limit 

XU  ■  Bode  plot  upper  phase  limit 

XL  ■  Bode  plot  lower  phase  limit 

XP  ”  Listing  switch 


Figure  10,  Format  Used  for  BI0DYN-80  TAPE 19  File  of  Transfer  Functions 
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BIODYN  T.F.  05/20/80 

SEATED  CREWMAN  -  VIEWING  DISTANCE 


Title 


15.13.57. 
-  9  CM 


D** 

RHD/DZP 


Dump  of  f ile 


DIODYN  T.E. 


BIODYN  T.F.  05/20/00 

SEATED  CREWMAN  -  VIEWING  DISTANCE 


15.13.57. 
9  CM 


Title 


RHD/DZP  "«■ - 

D** 

FIRST-ORDER  DIPOLES  CANCELLED 
46.384000 
100.00000 

SECOND -ORDER  DIPOLES  CANCELLED 
.44230900  .  8.3929300 

.60022000  f  50.002500 

.65024400  r  316.83100 


First  transfer 
function 

Cancelled  dipoles 
(only  exact  can¬ 
cellation 


numerator: 


-.81361 


(0 

» 

)  (0. 

)  (  10 

.955 

>  ( 

11.950 

< 

22.500 

)  (  107.75 

)  <  53 

5.79 

) 

<  < 

.22639 

.  a. 4561 

f 

1.9144 

f 

8.2365 

>>  ^ 

<  < 

. 78552 

»  13.624 

t 

10.702 

f 

8.4311 

>> 

(  ( 

.82761 

«  20.272 

f 

16.777 

T 

11.379 

>>  y 

<  ( 

-.16669 

f  21.651 

-3.6090 

f 

21.348 

>>  ( 

<  ( 

» 22759 

f  47.017 

* 

10.700 

f 

45.783 

>>  ) 

DENOMINATOR 

: 

( 

10.957 

)  (  11.983 

)  <  22 

.499 

)  ( 

107.81 

( 

531.02 

) 

(  / 

. 22869 

»  8.4602 

t 

1.9366 

r 

8.2437 

(  ( 

.40550 

r  13.701 

f 

5  •  55'j8 

r 

12.524 

(  ( 

.24026 

F  16.682 

r 

4.0081 

r 

16.193 

>>  ( 

(  < 

.82362 

f  20.305 

r 

16.723 

f 

11.516 

(  < 

•  94. >991*  - 

01.  31.293 

t 

2.9509 

f 

31.153 

>>  ) 

(  ( 

. 19299 

.  46.973 

r 

9.0600 

9 

46  091 

)>  ’ 

■ . 13057F- 

or.  - 

— 

>h- 


)\- 


High  frequency  gain 
First-order  zeros 

Second-order  zeros 


First-order  poles 
Second-order  poles 


Low  frequency  gain 


r:  1  u  I » !  N  I  .  F  . 


MODYN  f  ,F .  05/20/80  15.13.57. 

->l  AiCD  CREWMAN  -  VIEWING  DISTANCE  -  9  CM  *- 


RHLi/DZP 


XXX  MAGNITUDE  (DEO  - 

20.  0.  -20.  -40. 


-60. 


I -Hi. 


Vj. 


x  Magnitude 


-80. 


100. 


Title 

Top  scale 

XXX  is 
magnitude 


28. 


20. 


J.  4. 


I  0. 


I 


5.0 
4 . 0 

2  >>S 


i  .  1 


i  ■  0 

}•'/  s 


X  0 

ox 

0  X 

X 
X 

X 

X 


t 

Phase 

I 


o 

x  □ 

X  0 
X 

ox 

0  x 

0  X 
0  x 
0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

l) 

0 


+3608  phase  shift  for 
o  for  points  below  -l8o° 


Values  at  ">"  (in  rad/sec) 
and  at  — -  (in  Hertz)  are 
exact . 


1.00. 


0. 


X 

-100.  — £00.  -300.  -400. 

000  =  PHASE  <  DEG )  - 


Bottom  scale 
for  000  is 
phase 


b)  Printer -Genera ted  Frequency  Response  (Bode  Plot) 
Figure  11  (Continued) 


0 
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Title 


*•  •  w  •>  / 

'• 1  •••  / 


Listing  of  frequency 
(rad/sec),  amplitude  (dB) 
and  phase  (deg)  at  20 
evenly- spaced  increments/ 
decade.  (This  listing  is 
only  generated  if  XL  =  1.0 
in  the  CHOICES  file  for  this 
transfer  function) 


.  *  '-I  ■.  I 


3  3 t  /', 

;;;  o  .  i 

P.40  >  0 
;]  9  4  ...  0 


.i.  V  *  6 


..  /  I.  ♦  - 

;  7  2  v  j> 

i.  7:3  .  ? 
!.  7  -1 . 7 
i.  /":;.  i 


c)  Listing  of  Frequency  Response  Values 
Figure  11  (Concluded) 
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Note:  i>  denotes  a  typed  blank  space 


Output 

BIODYN  TF 

PIVIB 

TF 

Shoulder: 

JiXS 

Chi 

Inertial  horizontal  shoulder  deflection 

(S) 

BXP 

Chi 

S 

Horizontal  platform  deflection 

bxs 

m 

Chi 

Inertial  horizontal  shoulder  deflection 

BZP 

Ch2 

s 

Vertical  platform  deflection 

tzi 

tm 

Ch3 

Inertial  vertical  shoulder  deflection 

BXP 

Chi 

s 

Horizontal  platform  deflection 

JSZ1 

Ch3 

Inertial  vertical  shoulder  deflection 

BZP 

CH2 

s 

Vertical  platform  deflection 

Head: 

DTH 

a 

Ch3 

Inertial  head  pitch 

(H) 

BXP 

Chi 

H 

Horizontal  platform  deflection 

DTH 

m 

Ch3 

Inertial  head  pitch 

BZP 

Ch2 

H 

Vertical  platform  deflection 

DZH 

m 

Ch2l 

Inertial  vertical  head  deflection 

BXP 

Chi  1 

H 

Horizontal  platform  deflection 

DZH 

m 

Ch2 

m 

Inertial  vertical  head  deflection 

BZP 

Ch2 

L 

Vertical  platform  deflection 

Stick 

let 

Chi 

Longitudinal  stick  deflection 

(C) 

DXP 

Chi 

lc 

Horizontal  platform  deflection 

tet 

m 

Chi 

Longitudinal  stick  deflection 

DZP 

Ch2 

c 

Vertical  platform  deflection 

Assumptions:  Of  the  3  vibration  Inputs  possible,  2 
are  X  and  Z.  Assume  X  is  Channel  1,  Z  is  Channel  2 
for  the  above  table. 


Figure  12  BI0DYN-80  Transfer  Functions  Which  May  Be 
Used  as  PIVIB  Input 


1.  Specify  title  for  run 

2.  BDMOD  -  enter  the  biodynamic  response  module 

a.  FREQS  -  specify  all  frequencies  required  for 
blodynamic  analysis 

b.  VIBR  -  specify  vibration  input/source  transfer 
functions 

c.  BIOTR  -  specify  blodynamic  transfer  functions. 
This  Is  the  Interface  between  BI0DYN-80  and 
PIVIB. 

d.  BDOUT  -  specify  blodynamic  outputs  for  printing 

e.  BCOMP  -  given  all  information  specified  above, 
compute  and  print  biodynamic  quantities 

3.  PVMOD  -  enter  the  pilot/vehicle  module 

a.  Specify  various  P/V  parameters  in  state  vector 
form 

b.  VBINT  -  controls  communication  between  BDMOD  and 
PVMOD 

C.  PCOMP  -  perform  pilot/vehicle  computations  and 
print  results 

d.  FDREP  -  compute  and  print  frequency  domain  mea¬ 
sures,  e.g.,  describing  functions,  remnants, 
signal  spectra 

4.  Repeat  2  and  3  for  each  case 

5.  END 


B.  FOB  ASSISTANCE 


New  and  inexperienced  users  should  have  an  experienced  INTERCOM  or 
batch  user  assist  them  throught  the  first  application.  For  technical 
questions  concerning  this  program  or  models,  please  call  or  write: 
Henry  R.  Jex  at  Systems  Technology,  Inc.,  13766  S.  Hawthorne  Blvd., 
Hawthorne,  CA  90250.  Telephone  Number  (213)  679-2281. 
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APPENDIX  A 

DEFINITION  OF  PARAMETERS,  SYMBOLS,  RANGES  AND  TYPICAL  VALUES 

TABLE  1  PARAMETER  AND  SYMBOL  DEFINITIONS,  RANGES,  AND  TYPICAL  VALUES 


PARA¬ 

METER 

R3RTRAN 

wme 

UNITS 

LOCATION 

DEFINITION 

SOURCE* 

(Refs.) 

- REASONABLE - 

- 1 

TYPICAL  VALUES  (STD  CASES)  | 

RANGE 

MINIMUM 

MAXIMJM 

STDPLTt 

t 

STDTCRW 

LO£R  BODY 

* 

m 

kg 

Lower  body  ness  (Hips) 

T(l,2,3) 

10. 

25. 

14. 

14. 

Bb 

BB 

nAn/s 

Lower  body/seat  cush- 

F(l,4,5) 

100. 

2,000. 

1290.1 

315. 

ion  damper 

XKB 

NAn 

Lower  body/seat  cush- 

F(l,4,5) 

20,000. 

100,000. 

29585. 

49360. 

ion  stiffhess 

KfiS 

oa 

NAn 

Series  spring  g-adi- 

F(l,4,5) 

50,000. 

200,000. 

71519. 

119322. 

ent  in  lower  body 

0L 

TH. 

deg 

Seat  (hips)  tilt 

M(l) 

-20. 

+80. 

13. 

13. 

angle 

G 

G 

m/s2 

Trim  vertical  axel- 

0.1 

100. 

9.8 

9.8 

e-ation 

TORSO 

Mr 

XMT 

kg 

Torso  mass 

T(l,2,3) 

10. 

30. 

18. 

18. 

it 

ZIT 

kg.nC 

Torso  inertia 

T( 1,2,3) 

0.5 

2. 

0.8 

0.8 

&TB 

BFB 

N.m 

Torso  damping 

F(l,4) 

10. 

100. 

16.75 

40. 

kTB 

XKTB 

r®/5 

N.m 

Torso  stiffness 

G,F(1,4) 

100. 

2,000. 

500. 

500. 

* 

r3 

9t 

THT 

<feg 

Angle  of  XLT  fron 

M(  1,3,4) 

-30. 

+90. 

10. 

-15. 

vertical 

ILt 

XLT 

m 

Torso  c.g.  to  hips 

T,M(1,2,3) 

0.12 

0.40 

0.15 

0.15 

pivot  length 

Ptn 

THTN 

deg 

Angle  of  XLTN 

M,T( 1,2,3) 

-30. 

+90. 

5. 

-5. 

|ltn 

XLTN 

m 

Torso  c.g.  to  neck/ 

T,M(1,2,3) 

0.15 

0.30 

0.3 

0.3 

pivot 

P  s 

THS 

deg 

Angle  of  XLS 

M,T(1,3) 

-30. 

+90. 

5. 

-5. 

ILs 

XLS 

m 

Torso  c.g.  to  shoul- 

0.10 

0.30 

0.3 

0.3 

dsr  length 

f’STDPLl 

P*  =  STD  pil 

ot,  stiff  centerstick 

+"SDTC3W'  =  STD  craman,  hand  on  knees 


♦Soirees  Legend:  Measured  (e.g.,  photo)  or  specified  ty  situation  (continued) 

Tabulated,  cadaver  or  anthropometric  tables 
Fitted,  to  transfer  function  data 
TSjesstlmate  (on  basis  of  physical  properties) 
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TABLE  1.  (Continued) 


1 - 

1 

PAR fir- 

FORTRAN 

UNITS 

LOCATION 

SOURCE* 

RANGE 

TYPICAL  VALUES  (STD  CASES)  | 

METER 

(WE 

CEFINITICN 

(Refs.) 

MINIMJM 

MAXIM*!, 

STDPLT 

S1DCRW 

FEO( 

* 

m 

kg 

Neck  mass 

T(l,3,6) 

0. 

3.0 

0.0** 

1.24 

In 

m 

kg. 

Neck  inertia 

T(l,2,3,6) 

0. 

0.010 

0.0** 

0.0024 

Bn 

BNT 

N.m 

n37s 

N.m 

Neck/torso  danper 

G,F(1,5) 

0. 

0.40 

0.25 

0.1 

% 

XKNT 

Neck/torso  stiffness 

G(l,5) 

10. 

80. 

50.0 

50.0 

F3 

9n 

THN 

deg 

Angle  of  XLN 

M(l,3) 

-90. 

+90. 

-20. 

-30. 

Ln 

XLN 

m 

Neck  length 

M,T(1,2,3,6) 

0.05 

0.25 

0.1 

.12444 

lN1 

XLN1 

m 

Neck  c.g  to  neck/ 
torso  pivot 

T,GU,3) 

0.02 

0.15 

0.05 

.06222 

* 

m 

kg 

Head  mass 

T,G(J,2,3,5) 

2.0 

10.0 

helmeted 

4.34 

3.1 

IH 

ZIH 

kg  .m2 

Head  inertia 

T,G(1,2,3,6) 

0.010 

0.100 

helneted 

0.039 

0.0303 

Ml 

N.m 

rS7s 

N.m 

Head/neck  danper 

G,F(1,8) 

0. 

0.40 

0.0** 

0.126 

XKW 

Head/neck  stiffness 

G,F(1,4,8) 

10. 

30. 

10.** 

15. 

?a 

<h 

CH 

- 

Head/neck  conplianoe 

*1,8) 

0.0 

1. 

0.0** 

1. 

parameter  (CH  =  0. 
locks  head  on  neck) 

locked 

unlocked 

1 

TH 

deg 

Angle  of  H.G(p1vot-CG) 

M,T(1,2,3,6) 

^0. 

+150. 

70. 

60. 

ILt 

1 

l9v 

XLH 

m 

Head  c.g.  to  head/ 
neck  pivot  length 

T,M(1,2,3,6) 

0.0 

0.100 

0.0* 

0.025 

THV 

deg 

Angle  of  line  of 

Sl(Jlt 

M,T(1,3) 

-90. 

+90. 

-30. 

-30.  | 

1 

1 

IVd 

1 

VD 

m 

Viewing  distance 
from  head/neck 

HD 

0. 

1000. 

0.75 

0.6863  I 

1 

1 

1 

(pivot  | 

_ l 

*Sources  Legend:  **Head/Neck  Locked 

Measured  (eTg.,  photo)  or  specified  by  situation 
Tabulated,  cadaver  or  anthropometric  tables 
"Fitted,  to  transfer  firction  data 
TLesst  irate  (on  basis  of  ptysical  properties) 
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TABLE  1.  (Continued) 


TYP10L  WUK  (51D  CASES)  | 

PARA- 

FORTRAN 

UNITS 

LOCATION 

S01FCE* 

RANGE 

PtlfcK 

IW€ 

DEFINITION 

(Refs.) 

mm 

5TCTCT 

rsims 

ARMS  (Upper,  Lover) 

mi 

XM1 

kg 

Upper  arm  mass 

T,M(1,2,3) 

1.0 

4.0 

1.372 

1.372 

Hi 

n 

kg. m2 

Upper  arm  inertia 

T(l*2,3) 

0.01 

0.05 

0.012 

0.012 

k 

TI 

cfeg 

Upper  arm  angle 

M.T( 1,3,8) 

-30.0 

+90.0 

40. 

15. 

14 

D1 

m 

Upper  arm  length 

( »  4  ♦  4) 

M,T( 1,3,8) 

0.2 

0.5 

0.29 

0.29 

ILa/4 

R1 

fMn 

Upper  arm  c.g. 
(factional  distance 

T,G( 1,2,3) 

0.2 

0.6 

0.44 

0.44 

j* 

W 

kg 

Lover  arm  mass 

T,M(1,2,3) 

0.5 

3.0 

1.017 

1.017 

H2 

72 

kg-m2 

Lover  arm  inertia 

T,(l,2,3) 

0.01 

0.10 

0.0152 

0.0152 

\{k 

7E 

deg 

Elbow  angle 

M,T(1,3,8) 

60.0 

170.0 

145. 

115. 

14 

DP 

m 

Lover  arm  length 

(  =  4  +  4) 

M,T(1,3,8) 

0.2 

0.5 

0.305 

0.305 

14/4 

R2 

m/m 

Lover  arm  c.g. 
(fractional  distance) 

T,G(1,2,3) 

0.2 

0.6 

0.5 

0.5 

j9u 

GRIP  INTERFACE 

TIJ 

*g 

Grip  interface  angle 

W) 

-45.0 

+45.0 

0.0 

0.0 

IBi/Kl 

TIF 

s 

Grip  inteface  time 
constant 

G,F(1,8) 

0.0 

0.5 

0.01 

0.1 

lcI 

a 

m/N 

Grip  interface  com- 

H,G,F(1,8) 

0.0 

1000.0 

.0000558 

1.0 

1 

pliance  Ki-1 

fixed 

"free" 

♦Soirees  Legand:  Measired  (e.g.,  photo)  or  specified  by  situation 
Tabulated,  cadaver  or  anthropometric  tables 
Fitted,  to  transfer  f met  ion  data 
Ties st irate  (on  basis  of  physical  properties) 
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TABLE  1.  (Continued) 


— 

— 

1 

RANGE 

TYPICAL  VALUES  (STD  CASES)  | 

PARA- 

FORTRAN 

UNITS 

LOCATION 

SOURCE* 

1 

l€7ER 

tm 

CEFINITION 

(Refs.) 

MINIUM 

MAXIMJM 

STDPLT 

STICK 

— 

w& 

kg 

Stick  Aland  mass 
(referred  to  yip) 

M,F(1,8) 

0.0 

1.0 

1.31 

4.0 

Bs 

as 

N/m/s 

Stick  danper 

H,F(1,8) 

0.0 

5.0 

2. 

240.0 

*s 

XKS 

N/m 

Stick  cradient 

*1,8) 

0.01 

20,000.0 

13000. 

10000.0 

Qcs 

ocs 

- 

Stick  coipllanoe  para- 

M(l) 

0.0 

1.0 

1.0 

1.0 

neter  (CKSO  means 
that  XKS-’* 

"fined" 

•fined" 

Sx 

sx 

N/g 

Ebbuei^it  sensitivity 

to  Nx 

F,M,S(1,8) 

-5.0 

+10.0 

0.0 

0.0 

Sz 

sz 

N/g 

Bobvei^it  sensitivity 
to  Nz 

F,M,G(1,8) 

-10.0 

+10.0 

0.0 

0.0 

% 

THC 

(feg 

Stick  angle 

M(l) 

-20.0 

+120.0 

90. 

90. 

XLC 

m 

Stick  length  (pivot 
to  center  of  yip) 

M(l) 

0.0 

2.0 

0.61 

0.55 

*sc 

XKSC 

Units 

Stick  output  scale 

F,G(1,8) 

0.01 

20,000.0 

13000. (N/m) 

0.01 

m 

factor  (rescales  stick 
displacement  to  "spec- 

ified"  rnlts  (e.g., 
deg,  cm,  etc.) 

bAR 

BAR 

N/m/s 

Ann  rest  danper  (nor- 
mal) 

Arm  rest  stiffness 

H.G(1) 

0.0 

1.0 

o.ot 

t 

o.ot 

t 

kAR 

XKAR 

N/m 

M,G(1) 

0.0 

10,000.0 

o.ot 

O.Ot 

(normal) 

t 

t 

Bf 

BAT 

N/m/s 

Arm  rest  danpe- 
(tangential) 

G,M(1) 

0.0 

150.0 

0.0| 

o.oj 

*f 

XKAT 

N/m 

Ann  rest  stiffness 

G,M(1) 

0.0 

10,000.0 

o.ot 

o.ot 

(tangantlal) 

t 

T 

lBR 

XLEP 

m 

Elbow  to  am  rest 
distance 

M(l) 

0.0 

0.5 

0.0| 

O.OJ 

Atur 

ARM 

- 

Fraction  of  arm 

G(l) 

0.0 

1.0 

o.ot 

o.ot 

weiyt  on  arm  rest 

t 

♦Sarces  Legend:  Measired  (e.g.,  photo)  or  specified  by  situation 
Tabulated,  cadaver  or  anthropometric  tables 
fitted,  to  trarsfer  f met  Ion  data 
^estimate  (on  basis  of  physical  properties) 
tfrwest  disabled. 
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TABLE  1.  (Continued) 


PARA- 

l^TER 

FORTRAN 

NAAE 

UNITS 

LOCATION 

DEFINITION 

- 1 - 

SOURCE*  |  REASONABLE 
(Refs.)  |  RANGE 

TYPICAL  VALUES  (STD  CASES)  | 

NBJROKJSCULAR  SYSTEM: 

Actuation  Elements: 

1 

W 

_ 

Switch  for  central  or 

1 

M  |  0.0 

1.0 

0.0 

0.0 

force  activation  dy- 

1  arm 

head 

names  0  for  arm  con- 

trol  1  for  head/neck 

1 

control 

1 

Si 

SI 

- 

Overall  scale  factor 

F,G(1,8)  |  0.0 

200.0 

104.4 

50.0 

Kk 

XK1C 

WN 

Neu-onuscular  actua- 

F,G(1,8)  |  0.0 

1.0 

.00044 

.00 

tion  game 

1 

Tcc 

TCC 

s 

Neuromuscular  actua- 

F,G(1,8)  |  0.0 

0.2 

.0909 

.0909 

tion  game 

i 

cs 

TCS 

s 

Neuromuscular  actua- 

F,G(1,8)  |  0.0 

0.3 

0.089 

.089 

tion  lag 

1 

Ka 

XKAA 

N'N 

Gain  of  force  activ- 

1 

F,G(1,8)  |  0.0 

10.0 

1. 

1. 

ation  dynamos 

1 

^a 

ZA 

- 

Oarping  in  force 

F,G(1,8)  |  0.0 

1.0 

0.8 

0.8 

activation  dynamics 

1 

wa 

VA 

rad/s 

Natural  frequency  in 

F,G(1,8)  |  10.0 

20.0 

16. 

16. 

force  activation 

1 

dynamics 

T 

a 

TM 

s 

Time  delay  in  force 

F,G(1,8)  |  0.0 

0.10 

0.0 

0.0 

activation  dynamics 

1 

_ 

^Sources  Legend:  Meastred  (e.g.,  photo)  or  specified  by  situation 
“Tabulated,  cadaver  or  anthropometric  tables 
"Fitted,  to  transfer  function  data 
Tauesst  irate  (on  basis  of  physical  properties) 
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TABLE  1.  (Continued) 


PARA- 
ft  IER 

FORTRAN 

NWE 

UNITS 

LOCATION 

DEFINITION 

SOURCE* 

(Refs.) 

— mams. — 

RANGE 

mill 

mill 

STDPLT 

STDCFaf 

NBJROMJSCULAR  SYSTEM; 
Actuation  Elements: 
continued 

EH 

BM 

NAn/s 

"Hills  law"  darper 

G,F(1,7,8) 

0.0 

20.0 

1. 

1.0 

km 

XKM 

N.m 

Spring  in  naraius- 
cular  system 

G,F(1,7,8) 

0.1 

20.0 

2. 

2. 

Be 

I 

NAn/s 

Series  elastic  ele- 
nent  danper 

F(l,7,8) 

0.0 

10.0 

2.431 

2.43 

Kg 

XKE 

N/m 

Series  elastic  ele¬ 
ment  g-adient 

G,F(1,7,8) 

0.1 

100.0 

40. 

40. 

Bf 

BT 

NAn/s 

Tendon  darrper 

G,F(1,7,8) 

0.0 

10.0 

0. 

0. 

kt 

XKT 

N/m 

Tendon  g-adient 

6,F(1,7,8) 

0.1 

300.0 

80. 

80. 

Ksp 

XKS> 

N/m 

Hjscle  spindle 
model  gain 

G(l,7,8) 

1.0 

10.0 

5. 

5. 

Tsp 

TSP 

s 

Miscle  spindle  lead 
time  constant 

G(  1,7,8) 

0.02 

0.15 

1/11. 

1/11. 

T1* 

TPS 

s 

Mjscle  spindle  lag 
time  constant 

G(l,7,8) 

0.0 

0.1 

0.0 

0.0 

Tss 

TSS 

s 

Mjscle  spindle  hicfi 
frequency  lead  time 
constantly 

G(l,7,8) 

0.0 

0.1 

0.0 

0.0 

>XG 

N/N 

Golgi  tendon  organ 
model  gain 

G(l,7,8) 

' 

0.1 

1.0 

0.5 

— 

0.5 

T9 

IB 

s 

Golgi  tendon  organ 
time  constant 

G(l,7,8) 

0.02 

0.1 

1/18. 

1/18. 

h 

tz 

s 

Golgi  tendon  organ 
time  constant 

G(l,7,8) 

0.0 

0.1 

0.0 

0.0  | 

I 

TP 

TP 

s 

Golgi  tendon  organ 
lag  time  constant 

G(  1,7,8) 

0.0 

_ 

0.1 

0.0 

0.0  | 

1 

♦Sources  Legend:  Measired  (e.g.,  photo)  or  specified  by  situation 
Tabulated,  cadaver  or  anthropometric  tables 
"Fitted,  to  transfer  firction  data 
"Guesstimate  (on  basis  of  physical  prcperties) 
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TABLE  1.  (Concluded) 


i - 1  i - r 

— r 

- r 

— israwai — r 

- TYPICAL  VALUES 

~1 

iiii 

i 

i 

RAAGE  | 

(STAfOARD  CASES) 

1 

i  i  i  i 

I  PARA-  |  FORTH*  |  UNITS  | 

LOCATION 

i 

i 

r 

SOURCE*  | 

1  1 

1  1 

1 

1 

r 

i 

1  FCTER  I  N*C  |  | 

IIII 

Oaf  inition 

i 

i 

(Refs.)  | 

1 

MINIMUM  |  HVXIMJM  | 

1  1 

S1DPLT  |  SIDCRW 

1 

i 

i 

1 

1 

1 

1 

1 

1 

1 

1  IMAGE  FIXATION STIBJLO- 
1 - CCUL*  SERVO - 

1 

1 

I 

1 

1 

1 

1 

r 

i 

i 

i 

1 

1 

1 

r 

i 

i 

1  Koe 

! 

i 

i 

i 

| 

MODE 

1  rad 

1  rad 

1 

I 

1 

I  Fixation  error  gain 

1 

1 

1 

1  F  (9) 

1 

1 

1 

1 

1 

1 

1 

1 

1.0  | 

1 

1 

1 

10.0 

1 

1 

1 

1 

7.06 

I  7.06 

i 

i 

i 

i 

i 

j  *RE 

1 

i 

i 

i 

i 

XKRE 

1 

1 

I  rad 

1  rad/s 

1 

1 

I  Fixation  rate  gain 

1 

i 

1 

1  F  (9) 

1 

I 

1 

1 

1 

1 

1 

1 

0.0  | 

1 

1 

10.0 

1 

1 

1 

1 

1 

0.0 

I  0.0 

1 

|  i 

i 

i 

i 

I  TV 

1 

i 

I 

i 

TV 

1 

1  s 

1 

1 

i 

(  Tims  da  lay  in  fixation 

1  loop 

1 

!  F  (9) 

1 

1 

1 

1 

1 

0.0  | 

1 

1 

0.10 

1 

1 

1 

1 

0.045 

I  0.045 

1 

i 

i  i 

1  a 

I 

I 

i 

ALPHA 

1  rad/s 

1 

| 

1  "Pseudo  Integrator"  bre* 

1  frequency 

1 

1  G,  F(9) 

I 

| 

1 

1 

l 

0.0  1 

1 

1 

0.5 

1 

1 

1 

0.3 

1  0.3 

i 

i 

l  i 

|s 

l 

I 

i 

I 

XKP 

i  rad 

I  Fa? 

1 

1  T»-get  pirsuit  gain 

1  (Optional) 

1 

1  G  (9) 

1 

1 

1 

1 

1 

1 

0.0  | 

1 

1 

1.0 

1 

1 

1 

1 

0.0 

1  0.0 

1  j 

I  *IE 

r 

i 

i 

i 

*IE 

1 

1 

| 

I  Position  gain  frm 

1  "Vestibular"  sensor 
| 

r 

1  F  (9) 

1 

1 

— r~ 

i 

i 

l 

1 

0.1  | 

1 

1.0 

T 

1 

1 

i 

0.67 

1  0.67 

Jf 

j  i 

b 

i 

i 

| 

MOVE 

1 

1 

1 

1 

1  Velocity  gain  from 

I  "Vestibule"  sensor 

1 

|  F  (9) 

1 

1 

i 

l 

0.0  | 
j 

0.1 

1 

1 

I 

0.1 

j  01 

! 

i 

1  tc2 

i 

i 

JL 

TC2 

1 

1 

1 

i 

I  "Vestibular"  lag 

1  F,  T(9) 

1 

1 

i 

_i _ 

0.0051 

1 

0.015 

1 

1 

1 

0.01 

I  0.01 

i 

i 

i 

p 

i 

1  s 
| 

I  Ocala-  Serve  Lag  time 
constant 

fTT9] 

1 

i 

(ojrr 

1 

HIT 

T 

| 

in - 

I - O — 

| 

i 

1 

J  tlm 

i 

i 

i 

TLM 

1  s 

1 

1 

1  Ocular  Servo  Lead  tine 

1  constant 

1  F  (9) 

1 

1 

i 

i 

i 

0.0101 

1 

1 

0.020 

1 

1 

1 

1 

0.0125 

1  0.0125 

1 

1 

1 

1 

i 

i 

1  QE 

1 

1 

1 

1 

j 

1 

1 

1 

Switch  for  fast  mode  in 
Ocular  Servo 

1.  use  X,l£ 

0.  VE  =  infinity 

1  M 

1 

1 

1 

1 

I  0.0 

1  Pp”) 

1 

1 

1 

1 

1 

1 

1.0 

|  1.0 

1 

1 

1 

| 

1  1.0 

1 

1 

1 

| 

1  1 

1 

1 

1 

I 

1 

| 

jce 

1 

1  x 

| 

i  Ocular  Servo  Damping  Ratio 

1  F  (9) 

1 

I  0.2 

j 

1 

1 

1 

1.0 

|  0.65 
| 

|  0.65 

| 

1 

1 

j 

1  W 

1  O 

1 

1  V£ 

1 

I  rad/s 

|  Ocular  Servo  Natural 

1  F  (9) 

i  200.0 

1 

1 

400.0 

1316.23 

|  316.23 

| 

1 

I 

■  ® 

1 

1 

1 

1 

frequency 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

-1 

♦Sarces  Legend:  Iteasured  (e.g.,  photo)  or  specified  by  situ*1on 
T*ulated,  ca<hver  or  anthropometric  tables 
Titted,  to  transfer  function  data 
TLesstlwete  (on  basis  of  physical  properties) 
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TABLE  2. 

PARAMETER 

FILE 

STRUCTURE 

(SHORT  FORMAT) 

A. 

Lower  Body 

Mb 

Bb 

kb 

kbs 

eL 

G 

B. 

Torso 

Mp 

Btb 

®T 

% 

0TN 

%N 

0s 

Bs 

C. 

Neck 

% 

In 

BNT 

Knt 

®N 

% 

D. 

Head/Display  Viewing 

% 

% 

®HN 

kHN 

cH 

% 

% 

®V 

VD 

E. 

Arms  (Upper  then  lower) 

«1 

61 

h 

M2 

*2 

0E 

L2 

W1^ 

F. 

Grip  Interface 

eU 

bi/ki 

CI 

G. 

Stick 

Ms 

®s 

KS 

CKS 

Sx 

sz 

ec 

Lc 

Ksc 

43 


APPENDIX  B 

DEFINITION  OF  FORCING  FUNCTION  AND 
RESPONSE  VARIABLES 


TABLE  3  INPUT  (FORCING  FUNCTION)  VARIABLES 


1 

|  MATRIX 
|  COLUMN 
(  NUMBER 

1 

1 

INPUT 

VARIABLE 

NAME 

COLUMN 

CODE* 

UNITS 

POSITIVE 

DIRECTION 

DEFINITION  AND  USAGE 

1 

1 

1 

1  1 

1 

Nzp 

NZP 

g 

Upward 

Platform  acceleration 

1  2 

1 

**> 

DZP 

m 

Upward 

Platform  deflection 

1 

1  3 

1 

Nxp 

NXP 

g 

Forward 

Platform  acceleration 

1 

1  4 

1 

% 

DXP 

m 

Forward 

Platform  deflection 

1 

1  5 

1 

0FD 

N 

Forward 

Force  input  to  stick 

1 

1  6 

1 

Bzp 

BZP 

m 

Upward 

Special  Platform  deflection* 

1 

1  7 

1 

Bxp 

BXP 

m 

Forward 

Special  Platform  deflection* 

1 

1  8 

NMC 

NMC 

N 

Neuromuscular  command  within 

1 

CNS 

1  9 

1 

1 

9ti 

TTI 

rad 

Test  input  into  fixation  error 

♦Small  "b"  represents  a  typed-blank  space. 

Inertial  shoulder  (or  head)  motion 

♦These  are  used  to  get  transfer  functions  o7  Platform  motion 
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TABLE  4  OUTPUT  (RESPONSE)  VARIABLES 


MATRIX 

OOliMI 

NLMBER 

RESPONSE 

VARIABLE 

tm 

wvmix 

CCUJW 

CODE* 

UNITS 

POSITION 

DIRECTION 

DEFINITION  AND  USAGE 

1 

c 

fCfi 

m  or  N 

Forward 

Stick  output  (Kq  used  to  scale  as  foroe) 

2 

. 

fca 

FCA 

N 

Forward  on 
stick 

Interface  foroe  at  stick  yip 

3 

Aq 

HU 

ID 

Upward 

Vertical  shoulde*  deflection  relative  to 
platform  (zs-Zp) 

4 

A01  + 

on 

rad 

Pitch  ip 

Lowe-  ami  angle 

5 

A0a 

DTA 

rad 

Pitch  up 

Lower  amt  angle 

6 

fly 

F1Y 

rad 

Ip  an 
upper  anti 

Vertical  interface  force  at  elbow 

7 

N 

Muscle  nrodel  foroe  referred  to 
upper  am  c.g. 

8 

Li 

PLI 

N 

Forward 

Gr  interface  deflection  (stick  to 
lower  am) 

9 

fix 

FIX 

N 

Forward  on 
upper  arm 

Horizontal  interface  foroe  at  elbow 

ID 

*n 

|A)M 

m 

Internal  state  in  muscle  model 

11 

fi* 

FNS 

N 

Forward  on 
upper  am 

Horizontal  interface  foroe  at  shoulder/ 
upper  am 

12 

fvs 

FVS 

N 

Ip  on  upper 
arm 

Vertical  interface  foroe  at  shoulder/ 
upper  am 

13 

J6XS 

m 

Forward 

Horizontal  shouldr  deflection  relative 
to  platform  (xj-Xp) 

14 

a&t 

an 

rad 

Pitch  ip 

Torso  rotation 

15 

1 

frh 

FtH 

N 

Forward  cn 

Head/heck  interface  force  | 

1  16 

1 

A6fl 

CRH 

rad 

Pitch  ip 

1 

Head  rotation  1 

1  17 
| 

fvh 

FVH 

N 

Up  cn  head 

Head/neck  Interface  force  | 

1  18 

1 

fby 

FBY 

N 

Up  cn  Torso 

l 

Vertical  hips/torso  interface  force  I 

1 

1  19 

F* 

FBX 

N 

*  i 

Forward  on  |  Hips/torso  Interface  force  1 

1 _ 

hips 

1 

>  *No ter.  Small  y  represents  a  typed -bl arte  space. 
iNote:  Afy  (Column  4)  reans  small  signal  perturbation  cf  0\  given  In  fig.  3 
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TABLE  4  (Continued) 


MATRIX 

RESPONSE 

MATRIX 

POSITIVE 

COLUMN 

VARIABLE 

(DLtfN 

IANITS 

DIRECTION 

DEFINITION  AND  USAGE 

NUMBER 

NW€ 

OOOE* 

20 

* 

UEP 

m 

Up  the  seat 
beck 

Internal  state  in  hip  model 

21 

itp 

UBP 

m 

Up  the  seat 

Hip  deflection  relative  to 

back 

platform  (parallel  to  seat- 
back) 

22 

*1 

pa 

m 

Internal  state  In  rrusc  le 
nodal 

23 

a^N 

EfTN 

rad 

Pitoi  ip 

Neck  rotation 

24 

fnN 

FNN 

N 

Forward  on 

Neck/torso  horizontal  inter- 

neck 

face  force 

25 

^VN 

FVN 

N 

Up  on  neck 

Neck/torso  vertical  inter¬ 
face  force 

26 

Rrt) 

RH) 

tfeg 

Up  relative 

Head  point  of  Regard  relative 

to  displey 

to  a  display  on  the  platform  | 

27 

Ag 

(WA 

N 

fiiscle  rrodel  cormand  out  of  | 

spinal  cord  j 

28 

fa 

UFA 

N 

Corrmanded  force  out  of  force  | 
activation  dynamics  | 

29 

Faa 

FAA 

N 

Internal  state  in  force  activa-j 
tion  | 

30 

*9 

P>G 

N 

Output  of  gplgi  tendon  orgon  j 

sensors  | 

31 

* 

UPC 

N 

Spinal  cord  ccmrend  j 

32 

Ae 

P* 

N 

Error  between  neu’-onusculer  | 

ccmnand  and  proprioceptive  | 

stick  force  | 

33 

xsp 

XSP 

m 

Muscle  length  change  sensed  | 

by  spindle  nodel  | 

34 

Arn 

DRN 

m 

Relative  normal  deflection  of  | 
am  and  rest  I 

35 

Art 

CRT 

m 

Relative  tangential  deflec-  | 

tion  of  am  and  rest  | 

36 

AAi 

D3i 

m 

Vertical  conponent  of  head  | 

motion  relative  to  platform  | 

♦Note:  Small  *4"  r presents  a  taped -blank  space. 
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\jZ  . 


TABLE  4  (Concluded) 


MATRIX 

COLUm 

NUMBER 

RESPONSE 

VARIABLE 

NAME 

MATRIX 

COLUMN 

CODE* 

UNITS 

POSITIVE 

DIRECTION 

DEFINITION  AND  USAGE 

37 

M 

AC1 

rad 

Output  of  compensatory  action  in 
eye  fixation  loop 

38 

*E1 

AEI 

rad 

Output  of  "Pseudo  Integrator" 

39 

^MC 

TMC 

rad 

Comnand  to  Ocular  Servo 

40 

Ay 

bAV 

r/s 

Head  up 

Output  of  "Vestibular"  Sensor 

41 

eEH 

TEH 

rad 

Up 

Eye  rotation  relative  to  head 

42 

6M1 

TM1 

rad 

Internal  node  in  Ocular  Servo 

43 

e. 

THI 

rad 

Up 

Fixation  error 

44 

RED 

RED 

m 

Eye  up  rel. 
to  display 

Relative  eyem  point  of  regard  at 
display 

45 

eTC 

TTC 

rad 

Output  of  "Image  Pursuit"  opera¬ 
tion 

46 

ZT 

HZT 

m 

Up 

Display  inertial  displacement 
(vertical) 

47 

eE 

THE 

rad 

Up 

Inertial  eye  rotation 

48 

DXH 

m 

Forward 

Horizontal  component  of  head 
motion  relative  to  platform 

♦Note:  Small  "tf"  represents  a  typed-blank  space. 


APPENDIX  C 
EQUATIONS  OF  MOTION 

if  cKs  4  0 

(sin  ec)c  -  xs  -  L1(cose1)-^e1  +  L2(sinea)^6a  -  LjCcosejj)  =  - 

if  CKs  ^  0  if  C^g  =  0 

f - — - -  ^ 

|-(MgS2  +  BgS  +  Kg)  -  ^  -  M2(cosec)  s2  J  c/KgC  -  c 

+  (sinec)Lca  +  (cosec)lf)y  -  M2(cosec) (sineIJ)s2LI 
-  M2I^(cosec)  (cosea)s2^ea  +  (cosOcKB^s  +  KAR)^n) 


+  (cos0c)(Bfs  +  KF)(4rt) 

=  “  (  2/g  +  M2  cos0c)aZp  -  Sx/g  -  Fd 
3*  '  if^CKg  4  0 

( CO 80  ) 

Az1  +  L-,  (sino1  )/2i01  +  L^cosea)^  +  LjisinOjj)  +  — c  *  0 
k  |  I^2  +  (%M-,g  +  L-jMrjg  -  I'-,Favo)cos01  1^0!  -  I^j(sin0i  )fly 


+  Lb(cos01  )f^  +  La(cos01  )fns  +  La(sine1  )£vs  "  SiLaFM  -  0 


NOTE:  Favo  =  Arnip  (M1  +  M2)g  \  0a  =  90  +  0!  -  eEj 

=  eNZp  +  s2(^ep  +  BSp)  >  =  gN^t  +  s2(top  +  B^) 
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5 


|(I2  +  M2%cos20a)s2  +  (FaV()  -  M2g)Lesinea  -  I^fsinejf^ 


-  L2(cosea)f1y  -  La( sinea)fCa 
~  (Le  “  Ljb  +  Ldcos2ea)(BARs  + 

-  (l^sineacosea)(BFs  +  KF)£rT 

if  Cjr  4  0 

+  I^M2(cos0a) (cos0c)s2c  +  %(cos0a)M2(sin0jj)s2Lj  =  cosea)M2az 

p 

f1  +  fVs  ~  M1s2Az1  -  M1La(sine1)a2Ae1  =  M1az_ 

Jr  * 


CjCcosejj)^^  +  |cx(sin0Ij)aM2s2  +  Bj/Kj  +  1  ^  Lj 
+  CI(sin0IJ)(cos0a)M2Lds2^a 

-  CjCsine^f, 

- iW - 

+  CjCsiiiCjj)  [-  (Bfs  +  KF)(sinea)Art] 

=  CI(sin0IJ)M2aZp 

+  CjfsinOjj)  [-  (B^s  +  KAR)(cosea)Arn] 


if  Ck_  ^  0 

+  CI(sin0IJ) (cos0c)M2s2c 


9  if  cKs  J  0 

-M2(sinec)s2c  +  M2(coseIJ) s2^  -  |(M?Ldsin0a)s2  +  FaVQ  |  Aea 
Ksc 

+  fix  -  fCa  “  (BFs  +  KF)(cos0a)^rt  +  (B^s  +  K^)  (sinea)£rn  =  M2^ 
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fns  “  f1x  ~  MiS^s  -  M1La(c°s8i)s2^e1  = 

fvs  +  f VN  ~  fby  +  Mjfcose^s^p  -  MrLT(sineT)s2^6T  =  “  MraZp 

xs  +  (sineL)ubp  +  (Lpcosej  +  Lscos0s)A0t  =  B^ 

[LpS2  +  (B-pg  +  Bj^s  +  Kpg  +  Kppp  —  WE]  A0J  —  (Bj^ps  +  Kjjiji)A6jj 

—  ( Lscos9s  )fns  —  IjnjCcosajTjfJfjjj  —  Ls(sin0s)fvs 

-  LrN(sinaIiN)fYN  +  L]i(  cos9p  )fbx 

—  I»ji(sin9ji)fp)y  +  S^IiaFj^  =  0 

WE  =  ItC°s9t[(%  +  Mjj  +  Mip  +  M-|  +  M2)g  —  FavQ] 

+  Lscos0s[(Mt  +  Mg)g  “  Fav()]  +  LTNcosaTN(^I  +  *%)g 

^nh  +  ^%(s^n®L^s^  +  ^{(bpcosftji  +  Lpjjcos0rjijf  )s?A0kp 

+  MjjLj|(cos0jj)s2A9jj  +  MgLjp(cos0jp)s2A0g  =  Mga^ 

if  NM  0 

[(Ihs2  +  Bgjjs  -  %gLHcos0H)CH  +  KHN]^H_+_ 

-(BhnChs  +  kHn)A0N  ~  cHLH(sin0H)fvh  ~  CgLgCcosOi^fjjk  =  0 

fvh  —  Mjj(  cos0g)s2up)^  +  %(Lipsin0j  +  LjjySinQjijjJs^A&p 
+  %LN(sin%)s2A8N  +  MH%(sin9H)s^6H  »  ^zp 
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l8 


19 


(coseL)fby  +  (Bbs  +  %)ii^p  +  (sineL)fbx  +  MgS2^  = 

-  MbCcosSl)^ 

^bx  +  fnN  +  fns  -  MpI/rC  coso^s2/^,  _  ^(sine^Js2^^  =  -  MfaXp 


20 


/  BflS  +  Kb 

\  kBs 


U5P  ~  ’  bP 


0 


21 


(coseL)ubp  -  (Lrp3in0r  +  Lssin0g)A0T  -  Az-| 


23 

[%s2  +  Bjrf  +  KNT  -  +  )COS0N]A0N 

-  (sin©^)^  -  LN^coso^fjjj}  -  [(%  -  Ifli)sin0N  +  Ljjsin©^ 

-  [LjjCosBjj  +  (%  -  Iw^cosoj^fjjjj  +  [lHs2  “  J^S5Lbcos0jj]^0jj 

(BjjrpS  +  Kjqrji)^0rp  =  0 

24 

fnN  —  -^nh  +  *%(sin0L)s^ubp  +  Mjj(  I/jicos  9^  +  l/j^cos&jijj)s2Aftji 
+  %(LNlcos0N)s2A0N  =  %axp 


25 


fyjj  —  ^vh  —  f%(cos0L)s2ub  +  I%(Brsin9j  +  LjjjSin9pjj)s2A0rp 
+  %(I^1sin0N)s2A0N  =  %nZp 
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NEUROMUSCULAR  MODEL  EQUATIONS 


7 


Fm  +  (RjS  +  KT)La(Ae1)  -  (&J.S  +  Kj,)^  =  0 


10 


Pm  +  (Bes  +  K^Xm  -  (Bes  + 


0 


22 


Pm  -  (Bms  +  KM)La(A9T)  +  (Bms  +  %)xi  +  Pa  =  0 

if  NM  =  0 


if  NM  =  0 


27 


(Tpss  +  1  )(Aa  +  Ag  —  Ac)  ~  KSp(TSpS  +  1  )(Tsss  +  1 )xSp 


28 


29 


(Tas  +  2)FA  +  -  2)Faa  =  0 


F  —  F 
aa  raa 


This  set  to  1.0 
if  cua  =  0 


if  NM  =  0 


50  (Tps  +  1)A  +  Kg(Tgs  +  l)(Tzs  +  1 )%  =  0 


51 


(rcss  +  2)(Tccs  +  1  )AC  +  K1c(tcss  2)Ag  0 
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32 


NM  =  0 

/_ cos'  nc s‘’  |  c  —  sinecfca  +  (M2Iycoseacosecs2)AGa  -  cosoc  fiy 
\Ksc  ' 

o 

+  (MosinO|;jcosocs‘'')LI  -  (B^coseaCosepS  +  KARcoseacosec)'!Vn 
-  (BFsinoacosocs  +  KpsineacosPc)Art  +  Ap 
=  M2CosecaZp  +  NMC 
NM  «  1 

Ae  +  RHD/VT)  =  NMC 


33 


xsp  +  ~  xm  =  0 


26 


[I^sinfoT  —  Py)  +  I^jsin(0TN  —  9y)3^G,j  —  Vp/^GH  cos(0p  —  &y)u'5p 


+  Lfjsin( Ojj  -  Oy)^N  +  RHD  =  BZpcosOy  -  B^siney 


3^ 


A r  + 

^  n 


if  cKs  4  o 

c  cos(ec  -  oa) 
^sc 


+  (L2  -  +  Ljsin(ea  -  Gjj)  =  0 


35 


Art  -  c 


if  cKs  4  0 
^sc 


+  LjCOs( Ga  -  Bjj)  = 
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0 


36 


+  ( Iipsin0ip  + 


-  (cose^Ujjp  +  (%sineN)^eN  =  Bgp 

57  (2  +  V5^!  +  (Tvs  "  2^%Es  +  %e)6i  =  0 


38 

39 

40 

1+1 

42 

43 


(s  +  a)AEI  +  A  +  KxsAy  =  0 


eMC  "  0TC  “  AEI  +  %EAV  =  0 


(tC2s  +  1  )Ay  -  s^eH  =  0 
Qe  *  0  I  eEH  -  0M1  =  0 

%- 1  i  <'  *  s  *  S)e™ '  \  =  0 

e  ' 


(i  +  tems)8Mi  -  (i  +  tms)9Mc  =  o 


RED 

01  "Ip 


®TI ' 
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[LjSinfeT  -  9y)  +  Lj^sinCer^  -  9y)]^T  -  Vd^h 


— ubpcos(8L  -  &v)  +  ^eNLN8in(eN  -  ©v)  -  vD8eh  + 


Bz  cosev  “  B^siney 


RED 
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1*5 


Kp 

®TC  _  ZD 


0 


1*7 


eE  “  “  °EH 


0 


1*8 


Axjj  +  (sineL)ubp  +  (Lj,coseT  +  LrfjcoseTN)AeT  +  l^(coseN)^en 
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APPENDIX  D 


TYMSHARE  EXAMPLE  PROBLEM 


This  appendix  documents  the  use  of  CREATE,  BIODYN  and  PLOT  on  the 
Tymshare  System  31  PDP  10.*  The  computer  dialog  appears  on  the  next 
several  pages  for  a  typical  session.  An  existing  PARAMETER  file  is 
modified,  a  CHOICES  file  is  created,  BIODYN  generates  the  requested 
transfer  functions  (only  one  of  which  is  illustrated),  and  PLOT  produces 
a  quick-look  Bode  plot.  Throughout  the  dialog,  all  user  inputs/  respon¬ 
ses  are  underlined.  The  following  6  character  filenames  are  accessed  in 
the  course  of  this  example  (these  are  not  the  Standard  Pilot  and  Stan¬ 
dard  Crewman  which  are  available  on  Tymshare): 


File  Name  Title 


Existing  PARAMETER  file: 
Modified  PARAMETER  file: 
CHOICES  file: 

TF  file: 


BB19 

"Stiff  Stick 

SEMI 

"Semi-Supine 

CHOSAR 

TAPE19 

There  are  several  differences  between  the  CDC  version  of  the  BIODYN- 
80  package  and  its  Tymshare  counterpart.  These  are  listed  below: 

•  All  three  programs  run  Interactively  on  Tymshare;  all 
files  are  specified  by  the  user  during  execution. 

•  CHOICES  filename  can  be  an  existing  file,  whose  name 
is  to  be  modified.  (In  CDC,  the  CHOICES  filename 
must  be  a  new  name.) 


*A  user  planning  to  use  BI0DYN-80  on  Tymshare  should  be  familiar 
with  the  Tymshare  manual,  XEXEC,  especially  Section  3.  At  present, 
BIODYN-80  is  not  a  current  Tymshare  Library  Program,  and  a  potential 
user  should  contact  the  second  author  for  the  required  procedures. 


•  Because  of  the  above  differences,  there  are  slight 
changes  In  the  dialog  for  the  Tymshare  version  (e.g., 
the  NEW  PARAMETER  PILE?  query  is  not  used,  since  the 
PDP-10  software  can  determine  internally  whether  or 
not  a  new  filename  has  been  input). 
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D  •  1  Running  CREATE 


+ 1 


-LXL  LIGATE 


LEADING 
"Vf  CUT  ION 


Execute  CREATE 


INPUT  PARAMETER  FILENAME  J 

fci<19 


LISTING 


DESIRED* 


Interactive  input  of 
PARAMETER  file  to  be 
modified 


<LO>NG  OR  t'SH)ORT  LISTING7 

M 

PARAMETERS  F OR 2 
jTIFF  STICK 
LOWER  BODY 


I 1. 000  2290.1 

29595. 

?  1  5  l  ? . 

1 3 . 000 

TORSO 

18.000  0.80000 

J  6 . rso 

500 . 00 

10.000  0. I 3000 

2.0000 

0 . 30000 

5.0000 

NECK 

0.00000  0.00000 

-20.000  0.10000 

0.25000 
0.50000C  01 

50 . 000 

HEAD/ DISPLAY  VIEWING 

1.3400  0.39000F *01 

0.00000 

l  .  0000 

0  .C-'OOO 

70.000  0.00000 

-30 . 000 

0 . 75000 

ARM  (UR'  FRf  LUWLK  > 

1.3720  0.,  12000E  -01 

•10.000 

0 . 290 AO 

0  11000 

1.017'  0*.  15/OOF  01 

i  v*.  .oo 

0 . 30500 

0.50000 

DR IF  INTERFACE 

0.0000  -  O.K'OOOt -01 

0.2573  IE -01 

\>r  icf 

0.31000  2.0)0 A 

90.00“  0. i 1000 

13900. 

13900. 

l  .0000 

0,00000 

ARM  F\ !”  5  T 

O.OOA00  O.'MOOO 

0 . 00000 

0.00000 

0,00000 

/ urgmi  - r  ijl  ap  .  r  s  r km 

o  01. 40 

o.;:oi3*r-oi  o.vovo?e-oi 

0.3938RE  -01 

1. 0000  •S.t'0000 

16.000 

0 . 00000 

l  .oo;..-)  -.  O00 

7.4312 

40.000 

0.00000 

D.oooo  o.'.  ovoof-oi 

A . OAOOO 

0 . 00000 

0,50 >00  A.55554F-01 

0,00000 

0.00000 

image  ;  rxATTnrj/vesriRin 

i.l-IJi  III  Ak  GtRVl) 

7. 0.*vr.  O.'-OOOO 

*) ,  A 7000  0.  10000 

0.  I5000E-01 
0. 1 AOOOE  0 1 

>.  10000 

•'.00000 

o.ioo  ■;  :  .oot -<  i 

1  .yOOO 

•  • . 65000 

<16.73 

CHANGES  DESIRED  - 

YC 


INMtl  MMFMONJf  •  5  CHA*5> 
r-j  -  xt  r ,  t  yf  r  yxx . 


no. 


AND  VAMif  FOR  f  AC!H  CHANGE  t 


9.8000 

0.30000 


l  Generate  short  listing 

f  of  file  to  modify 


0.00000 

0.00000 


30.000 


/ 

l —  Make  changes  to  file 


(continued  on  next  page) 
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D.1  (Continued) 


■  Ui.'M.  I.rslKM-' 


a  ijk  I’jH-iihr  i.  ;sriri!',? 

ri 

■  . .  M  t-  II  I :  ■ .  f  TIST 
;  ill  -  .ill-  INI. 


j  .i  r  pin 


:  .  ■"■o 

iT'-'O.  i 

-*5r.!3 . 

’IE  19. 

G  r  ■ '  n 

U  .  -->00 

o.neooo 

I  *  .  *’V. 

'•OO.vO 

. 0  DO 

r: . 

. .0000 

0 .  .50000 

. n ' '00 

■:* .  .loo 

.  .1,00  0 

•  0.  wOO 

■>  j 

0.  IC'Vm. 

C  OOO  ■•!:  -01 

E  r M-r 

1  .V 

0 !  1  u:  i  10 

; .  : ; ;  o 

0  •  ‘  90-»h  -0  1 

I . OoOOO 

1 .0000 

c .  :.oooo 

no.  oo 

0 .  C/'i./O0 

: 0.000 

0.  7C.000 

'■  id  I..T  FI- 

r  . 

I  '.lUl-  i'. 

•  .  ’'’.'O 

.  i 'OCC-  01 

’•.000 

0 . J9000 

0. 1 1--0 

1  .  /  1  "  ’• 

.  i  ..  m 

:  r.,.oo 

0.  50r,(:0 

•  <>o 

r  i  r-  :  :i  r 

.  t  ■ 

r* 

.  ■-wOO 

■.10-  1  • ! 

o.o  .  'fur  0  1 

.  I'OO 

I  .*900. 

1 .0000 

0 ♦ 00000 

0 .  M  ’  •- 

:  7v--o. 

r  !  -iT 

.  :r  i  v ; 

in 

»•  ’  ,  :  1  1  M 

•  OOOOO 

0.00000 

0 . 00000 

»  i . 

10 

.  ;  <  *r 

1 

o.  -01 

O.HVHBldi:  01 

} . 

C 

10.000 

v).  00000 

!  o:;oc 

’.0 

• . 

10 . ooo 

0 . 00000 

"• .  0000 

I'.VO  "-Of  -01 

•  .00000 

0.00000 

.  ,  >'>GC 

/  M  -01 

.00000 

0.00000 

i  « . 

i'i  \ 

,,i/o  r 

0  I.UJLY.R  :;EI-:00 

'  .  OOf.  GO 

o.  YJOiK-F  01 

0. JOOOO 

0 . 00 ooo 

. <>•>: 

0.  ('/•/('  > 

'. 1 0000 K -01 

0.  1  "  "'1  01 

l .0000 

•>.  ATiOOO 

51ft.  JJ5 

r :  : 

'thf 

fOk  r.iKOhK 

n  a  f  n  e  :  v 

Modify  title  for 
PARAMETER  file 


Short  listing  of 
modified  file 


Interactive  input 
of  filename  on  which 
to  store  modified 
PARAMETER  file 


(Continued  on  next  page) 
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D.1  (Concluded) 


INPUT  CHOICES  FILENAME 
CHOSAR 


BIODYN  TFS  DESIRED  FDR  PIVIB 


M. 


’h 


transfer  function  input  : 

FIRST  LINE-RESPONSE  MNEMONIC  *  FORCING  FUNCTION  MNEMONIC 
<  AAA  * AAA )  *  ENTER  XXX  TO  STOP 
SECOND  LINE-PLOTTING  INFORMATION*  5  ITEMS  : 

DOPE  LOWER  FREQ.  LIMIT 
BODE  UPPER  FREQ  *  LIMIT 
BODE  UPPER  PHASE  LIMIT  <0. 

BODE  LOWER  PHASE  LIMIT  (0. 

LIST  <1.  TO  LIST  TABLE*  0. 

IF  NO  PLOT  DESIRED*  fcNTER  0. 


DEFAULTS  TO  200.) 
DEFAULTS  TO  -400.) 
FOR  NO  LIST) 

FOR  ALL  ITEMS 


1. *100, *0. *0. *0. 
l.»100.»0.»0.»0. 

JSEEIE 

i..ioo.»o.»o.»i. 


MX. 


Interactive  input  of 
CHOICES  filename 

No  PIVTB  TF's  desired 


Format  for  CHOICES  file 
information  input 


Input  CHOICES  file 
information 


EXIT 


Exit  CREATE 


D.2  Running  BIODYN 


1 


EXE  PDN80fN7RFN 

nmm - 

EXECUTION 


Execute  BIODYN 


ENtCR  INPUT  FILE  NAME :  SEMI 

ENICR  CHOICES  FILE  NAMES  CHQSAR 

1 1  -Nar-80  ll!53 


CASE  S  SEMI  -SUPINE 

ENTER  SYSTEM  OUTPUTS  FILE  NAME'-*  TAPE 20 
NFU  TILE 

F-  DAT  I  NO  UNDERFIOU  PC  - 1 1  1246 

FLOATING  UNDERFLOW  PC- 111? 44 


Interactive  input  of 
PARAMETER  and  CHOICES 
files 

BIODYN  execution 
Case  title 

Interactive  input  of 
TF  file 


denominator: 


104446-26 

( 

16.591  > 

<  46.384 

( 

100.00  ) 

(  285.86 

( ( 

. 44231 

f  8.3929 

<  < 

.23697 

f  11.309 

<  ( 

.20726 

f  13.642 

( < 

. 19865 

f  24.860 

( l 

.32015 

f  32.898 

<  ( 

.83789 

f  32.931 

<  < 

.19448 

f  57.685 

<  < 

. 65024 

f  316.83 

. 1 8027E+1 1 


)  (  56.079  )  (  97.512 

)  (  333.22  ) 


3.7123 

,  7.5273 

)  ) 

2.9060 

f  10.929 

)  ) 

2.8274 

f  13.346 

)  ) 

4.9386 

f  24.365 

>  > 

10.532 

f  31.167 

)  ) 

28.251 

.  16.921 

)  > 

11.219 

f  56.583 

)  ) 

206.02 

f  240,70 

)  ) 

NUMERATOR J  DTH/DZP 


, 29296E* 27 

«  .28049 

<  . OOOOO  ) 

<  .00000 

>  <-10. 

621  > 

16.590 

(  46.384  ) 

<  97.067 

)  <100 

.00  ) 

-138.57 

<  209.05  ) 

<(  .44231 

f  8.3929 

3.7123 

f  7.5273 

>  ) 

(<  .21687 

f  11.488 

2 . 49 J  3 

,  13.714 

)  ) 

((  .98984 

.  16..J17 

16>»46 

f  7.39  16 

)  t 

' <  . US 1 30 

f  33.163 

28.232 

,  1/.100 

)  ) 

( (  . 20646  ♦ 

f  53.410 

11 . 440 

,  34.216 

)  ) 

<(  .63024 

f  316.83 

206.02 

f  740.70 

>  ) 

\  (  .  V6  4  5 

*  3Ja . 66 

.325 . 4 1 

.  84.321 

) ) 

.67996FF08 

.37719E 

•02 

mu  -123.31 

U-rtir-80 

1 1  .*  58 

y 


> 

> 


EXIT 


Printout  of  DTH/DZP 
transfer  function 


(See  Figure  9  for 
explanation  of 
format ) 


Exit  BIODYN 
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D.3  Running  PLOT 


fc  xf  rTFLQT>jCCUfc^C&OPXL 
;.'0At»  I  NG 
F*E'i:iiriON 

ENTER  TRFN  SYSTEM  TILE  NAME  J  TAPE20 

rl  •  f  MR  11 -Mar-80  11153 

STMT  -  Slip  1  HE 

I'M 

c:  /NZf' 

M/DZP 

riTH/DZP 


Execute  PLOT 

Interactive  entry 
of  TF  file 
Case  title 

Dump  of  TF  file 
contents 


TOPE  jo 


PLOT  execution 


,i .  rhP  II  Her  -  no  1  i  *c*^ 

•  Ml  -SUP  INK 

UTM/DZP 
l»*  *• 

r  ll:ST  tiM<ER  [il’  Oll  S  CANCEL  I.LU  » 

14 .  J-V  >000 
1 00*00000 

»F.rOMI’- ORDER  l'IPOLl.5  CANCfcLLED  i 
.  VI.’ ’••'WOO  *  il  •  .59 2*'  300 

.  1.Z fOJttOO  *  516.83100 

Nl/MERA  I  HI  ' : 


Transfer  function  DTH/DZP 

Identical  roots  in  numera¬ 
tor  and  denominator  are 
cancelled}  close  pairs  are 
retained 


.  rno-»v 

,0  5oO\>  ) 

4  . 00000 

v  7  > 

< -138.57 

(  . .?  1  43  7 

[  1  .  483 

<(  .98784 

t'. 817 

<  .  i’52  JO 

3. 143 

<  . J044A 

55.410 

(•  ,'Ma'j7 

334 . A A 

ri[  hum  i  tic.  1  nr. : 

•:  t-s.'.vt 

<  35  5  • .  ‘  > 

.  06.0/9 

•  t  .  - :  -  7 

.  11. 50V 

;<  .:o/:a 

,  13. 6 43 

<  •;  .  !V.-'A*..» 

.  71,860 

<(  .  ;.roi  5 

»  3J.898 

<<  .  h'j  ;nv 

.  52.731 

< <  .1/148 

.  5  7 . 48  '* 

,  •  VlVt  0.' 


>  <10.621 

)  < 

16.590 

)  <  ’09.05 

) 

2.4915  * 

11.214 

>  ) 

16.616 

2.391? 

) ) 

28.237  » 

17. 100 

)  ) 

11.110  ' 

51.216 

)  ) 

325.11  » 

26.321 

)  ) 

>  <  V  7 

517 

)  < 

285 . 86 

2 . 9040 

10.929 

>  ) 

2.8274 

13.544 

>  ) 

1.9486 

t 

74.365 

1  ) 

10.537 

t 

31 . 167 

>  ) 

20 . 25 1 

1 6 • 971 

)  ) 

11.21' 

56 .503 

)  ) 

> 


) 


Transfer  function  to 
be  plotted  (see 
Figure  1 1  for  explana 
tion  of  format) 


/ 


(Continued  on  next  page) 
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'*»**■■  • . 


D.3  (Continued) 


M.THP  U-Mrfr-80  11553 

SEMI  SUPINE 


DTH/D2P 

Vtt 


XXX  =  MAGNITUDE  <PB) 


40.  20. 


0.  -20.  -40.  -60.  -80. 


28 .  >  X 

X 


X  0 

X  0 

X  0 

X  0 

X  0 
0  X 

0  X 

0  X 


7.1  & 


5.0  0 

0 

4.0  0 

0 

0.5 -  0 

2.8  0 

0 
0 

2.0  U 

n 

u 

1.4  0 

0.2 -  0 

O 

1.0  0 


200.  LOO. 


0.  -100.  200.  -300,  100. 


000  PHASE  CDI.(i> 


(Continued  on  next  page) 


Bode  plot  of  DTH/DZP 
XXX  is  the  magnitude 


Frequencies  which  are 
denoted  as  hertz  (HZ) 
or  radians /second 
(R/s)  are  exact  at 


respectively 


000  is  the  phase 
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D.3  (Concluded) 


T  AF'K  70 

m.imf  n  Mur- flo  11:53 

:-T  hi -  SUPINE 


PTM/DZP 

r»$* 


h'/S 

I. 00 

I  .12 
3  .  26 

1 .41 
l  .53 
1.  "8 
2.00 
2.24 
2 . 5 1 
2.82 
3.1* 
3.55 
3 .  98 
4.47 
5.01 

5.42 
4. 31 
7.08 
f  .  94 
3.91 

10.00 

II. 22 
12.59 
14.13 
i  5 . 85 
l  7.  73 
1 9  ♦  y  5 
.•7.39 

70 .  I  3 

1  i  .  42 
55.4  8 
*9.31 
4  4.*7 
■ .  r. ,  *2 
',*.73 
v'.  10 
79 
»V,  43 
09.  13 
« 00 . 00 
xi  r 


r»* 

-48.34 

46.31 
44.27 

-42.21 

-40.15 

-38.07 

-35.96 

-33.83 

-31.67 

-29.46 

-27.20 

-74.88 

-22.48 

-19.97 

-17.34 

-14.56 

-11.58 

-8.38 

-4.91 

-1.17 

2.85 

7.30 

12.34 
16.45 
18.74*' 
30.72 

23.31 

26.34 

23.34 

77.92 

26.04 
23.40 
70.28 
17.12 
:  4.43 
I  2 . 60 
!  >.32 
3.59 
6.37 
1 . 35 
2.51 


ICG 

175.9 

175.4 

174.9 

174.3 

173.5 

172.8 

171 .9 

170.8 

169.7 

168.4 
167.0 

1 65 . 3 

163.4 
161.2 

158.6 

1 55 . 5 

151.7 

147.1 

141.1 

133.2 

122.9 
ItO.O 

89.2 

56.6 

24.4 

-J  .7 
-27.3 
-60  •  1 
-104.6 
-150.8 
-189.7 
-221.9 
-246.8 
-264.3 

275.8 
-  205 . 8 

299,4 
-31  1  .4 
-320.3 
-327.3 
-333.2 


Listing  of  frequency  (rad/sec), 
amplitude  (dB)  and  phase  (deg) 
at  20  evenly-spaced  increments 
per  decade  for  DTH/PZP 

(This  listing  is  only  generated 
if  XL  =  1 .0  in  the  CHOICES  file 
for  this  transfer  function) 


Exit  PLOT 
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APPENDIX  E 


INTERCOM  EXAMPLE  PROBLEM 


This  appendix  presents  a  typical  terminal  session  in  the  Intercom 
4.7  operating  system,  accessing  the  CSA  mainframe  at  Wright-Pat terson 
AFB's  ASD  Computer  Center.  The  dialog  is  annotated  so  that  a  typical 
user  will  easily  understand  the  basic  sequences  of  parameter  entry  and 
Job  steps.  All  user  responses  are  underlined;  each  is  terminated  by  a 
carriage  return. 

This  particular  session  was  an  exercise  using  a  9  cm  viewing  dis¬ 
tance  and  the  standard  crewman  (STDCRW)  parameter  set.  It  was  one  of  a 
series  of  runs  which  attempted  to  determine  optimal  display  distance  for 
minimizing  image  motion  of  a  vertically  vibrating  crewman.  The  example 
is  carried  far  enough  for  the  potential  user  to  see  how  the  programs 
Interface  with  INTERCOM  job  control  and  file  management  commands. 

The  steps  followed  in  the  investigation  are  listed  below: 

1)  Log  in  to  Intercom. 

2)  Attach  the  PARAMETER  file  to  be  modified  and  name 
it  TAPE20. 

3)  Attach  the  CREATE  program,  called  EXECRT . 

4)  Run  EXECRT,  make  changes  to  existing  file,  assem¬ 
ble  new  CHOICES  file. 

5)  Assemble  batch  job  to  run  BIODYN  and  PLOT,  using 
the  two  PARAMETER  file  (old  and  modified)  and  the 
new  CHOICES  file. 

6)  Submit  batch  job  to  input  queue. 

7)  When  job  completed,  list  the  output  file. 

The  user  is  advised  to  retain  in  his  permanent  files  only  those  which  he 
wishes  to  use  in  the  future,  and  to  frequently  purge  his  directory  of 
unneeded  parameter  and  output  files  in  order  to  forestall  job  failures 
resulting  from  file  space  overload. 
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Typical  INTERCOM  Terminal  Session 
user  responses  underlined 


ASP  COMPUTER  CENTER  INTERCOM  5.0 
SYSTEM  CSA 

'  PATE  05/20/00  T  T ME  1 4 . 4/, .  3V  . 

PLEASE  LOGIN 

ENTER  3-DIGIT  TERMINAL  T Ti  — 

05/20/80  LOGGED  IN  AI  14.44.52. 
WITH  USER  IP  CK 
EQUIP/PORT  12/001 

COMMAND-  REQUEST, TAPE7.HPF 

COMMAND-  REQUEST. TAPES. *PE 


1 .  Log  in  to  INTERCOM 


COMMAND-  ATTACH ,  T  AEE20 .  STOCKUJ 
PF  CYCLE  NO.  -  001 
COMMAND-  ATTACH.EXECRT 


PFN  IS 
EXECRT 

PF  CYCLE  NO.  -  001 
COMMAND-  EXECRT 


b 


2.  Attach  files  to  I/O  units 


3.  Attach  EXECRT( CREATE)  program  file 


NO 

YE 

LI 

LO 

YE 

NO 


L.  Execute  EXECRT  program  module 
new  file  v  j _  5*  Select  existing  data  file  (STDCHW) 

LISTING  DESIRED?^.  List  data  file 
(LO)NG  SR  (SH)ORT  I  TPTfNA'.'j 
ILOLNG  OR  (SH)ORT  '  tsriNl.vl 


7*  Program  rejects  unrecognizable 


NEW  CHOICES  FILE 


V 


BIODYN  rrs  DESIRED  FOR  PIVI 


“  "h 


8.  Long  listing  option  selected 

9.  Affirmative  on  CHOICES  option 

10.  Negative  on  PIVIB  option 


TRANSFER  FUNCTION  INPUT  ! 

FIRST  LINE-RESFONSE  MNEMONIC,  FORCING  FUNCTION  MNEMONIC 
(AAA, AAA)  i  ENTER  XXX  TO  STOP 
SECOND  LINE-PI. OTTTNG  INEORMA I  ION ,  5  ITEMS  1 
BODE  LOWER  FRECI.  LIMIT 

POPE  UPPER  FREQ.  LIMIT 

BODE  UPPER  PHASE  LIMIT  <0.  DEFAULTS  TO  200.) 

BODE  LOWER  PHASE  LIMIT  (0.  DEFAULTS  10  -400.) 

LIST  (1.  TO  LIST  TABLE,  0,  FOR  NO  LIST) 

IF  NO  PLOT  DESIRED,  ENTER  0.  POP  All.  ITEMS 

~FTNP  not  permissible,  pi  easf  rf INPUT }- 1 1  •  Illegal  parameter  selected 
RHD.PZP  1 

.  1 . 200  ,  ,  0  .  ,  0  .  ,  1  « 


MAX  FREQUENCY  RANGE  IS  .1  DECADE:' 
PLEASE  .m'  INPUT  f  NT  IRE  LINE 


12.  Program  screens  for  out-of-limit 
values 


1 1 1  99 •  r.9 1  ’ 0 1  ’ 1  - - -  13.  Finally,  valid  input 

■jLU - ’it.  Exit  from  CREATE 

STOP 


CP  SF  COWlTS  execution 


Session  continues  with  batch  run 
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Sample  INTERCOM  Batch  run 


COMMAND-  EDITOR 


H-  i 


.  . CREATE. S 


ENTER  LINES 

SAR »  CU15000C  ,  S  fCSA « 


L8Q00 15. RIEDEL. (213)679-2231 


ATTACH. TAPE? . CHOICE 1 . 

ATTACH, TAPF3.UIFU9, 

ArfACH.EXEBIO. 

ATTACH . EXEPLT • 

EXEB 10 ( TAPE A > OU  TPUT . TAPE 7 . TAPE8 > TAPE19. TAPE21 ) . 


'fcfcuiND.fSPEIT; 

EXEPLT (OUTPUT, TAPE19) . 
«EOR 


2 


.  .SAVE.GQFILF .NOSEQ 


•  * STORE , GOP IL£ .LU00015 

CT  ID*  L300015  PFN-GOF ILE 
CT  CY  =  001  0000012S  WORDS.: 

•  .end 


COMMAND-  BATCH, uOFILE. INPUT, HERE 


3 


h 


COMMAND-  FILES 


—LOCAL  FILES  — 

•VEXECRT  *00:-'ILE  *  INPUT 

HTAPE3  HTAPE20 
—REMOTE  EXECUTING  JOBS  — 
SARCK02 
COMMAND- 

TASK  , TEL  *  TERMINATED 
TASK  EXIT  WITH  OUTSTANDING  10 
PC  =  10664.4 
F3- 170000 
R0-:003770 


$6UTPUT  *TAPE7 


Note:  Output  files  generated  by  BIQDYN,  vnen  listed, 
are  identical  in  format  to  printouts  shown  in 
Appendix  D. 


.  Invoke  INTERCOM  editor 


.  Create  batch 
job  stream 


.  Preserve  file  containing 
Job  stream 

.  Batch  job  stream  to 
input  queue 


Job  terminated  by  system 
crash 
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